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Relativistic hydrodynamics has been extensively applied to high energy 
heavy-ion collisions. We review hydrodynamic calculations for Au+Au 
collisions at RHIC energies and provide a comprehensive comparison 
between the model and experimental data. The model provides a very 
good description of all measured momentum distributions in central and 
semiperipheral Au+Au collisions, including the momentum anisotropies 
(elliptic flow) and systematic dependencies on the hadron rest masses up 
to transverse momenta of about 1.5-2 GeV/c. This provides impressive 
evidence that the bulk of the fireball matter shows efficient thermal- 
ization and behaves hydrodynamically. At higher px the hydrodynamic 
model begins to gradually break down, following an interesting pattern 
which we discuss. The elliptic flow anisotropy is shown to develop early 
in the collision and to provide important information about the early ex- 
pansion stage, pointing to the formation of a highly equilibrated quark- 
giuon plasma at energy densities well above the deconfinement threshold. 
Two-particle momentum correlations provide information about the spa- 
tial structure of the fireball (size, deformation, flow) at the end of the 
collision. Hydrodynamic calculations of the two-particle correlation func- 
tions do not describe the data very well. Possible origins of the discrep- 
ancies are discussed but not fully resolved, and further measurements to 
help clarify this situation are suggested. 
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Hydrodynamic description of ultrarelativistic heavy-ion collisions 3 

1. Introduction 

The idea of exploiting the laws of ideal hydrodynamics to describe the 
expansion of the strongly interacting matter that is formed in high en- 
ergy hadronic collisions was first formulated by Landau in 1953P Because 
of their conceptual beauty and simplicity, models based on hydrodynamic 
principles have been applied to calculate a large number of observables for 
various colliding systems and over a broad range of beam energies. How- 
ever, it is by no means clear that the highly excited, but still small systems 
produced in those violent collisions satisfy the criteria justifying a dynam- 
ical treatment in terms of a macroscopic theory whic h follows idealized 
laws (see Section ETT|) . Only recently, with first datd^l from the Relativis- 
tic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory came 
striking evidence for a strong collective expansion which is, for the first 
time, in quantitative agreement with hydrodynamic predictionspl both in 
central and non-central collisions (in this case Au+Au collisions with center 
of mass (cm.) energies of 130 and 200 GeV per nucleon-pair) . 

The validity of ideal hydrodynamics requires local relaxation times to- 
wards thermal equilibrium which are much shorter than any macroscopic 
dynamical time scale. The significance and importance of rapid thermal- 
ization of the created fireball matter cannot be over-stressed: Only if the 
system is close to local thermal equilibrium, its thermodynamic properties, 
such as its pressure, entropy density and temperature, are well defined. And 
only under these conditions can we pursue to study the equation of state of 
strongly interacting matter at high temperatures. This is particularly inter- 
esting in the light of the expected phase transition of strongly interacting 
matter which, at a critical energy density of about 1 GeV/fm 3 , undergoes 
a transition from a hadron resonance gas to a hot and dense plasm a of 
color deconfined quarks and gluons. Lattice QCD calculations indicate^*' 
that this transition takes place rather rapidly at a critical temperature T cr j t 
somewhere between 155 and 175 MeV. 

In this article we review and discuss data and calculations which pro- 
vide strong evidence that the created fireball matter reaches temperatures 
above 2T cr it and which indicate a thermalization time below 1 fm/c. Due 
to the already existing extensive literature on relativistic hydrodynamics, 
in particular in the context of nuclear collisions, we will be rather brief 
on its theoretical foundations, referring instead to th e availab le excellent 
introductory materia l^ l ^ l and comprehensive reviews 

In Section [21 we briefly present the formulation of the hydrodynamic 
framework. Starting with the microscopic prerequisites, we perform the 
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4 Hydrodynamic description of ultrarelativistic heavy-ion collisions 

transition to macroscopic thermodynamic fields and formulate the hydro- 
dynamic equations of motion. These are equivalent to the local conserva- 
tion of energy and momentum of a relativistic fluid, together with continu- 
ity equations for conserved charges in the system. The set of equations is 
closed by providing a nuclear equation of state whose parameterization is 
also presented in that section. This is followed by a discussion of the initial 
conditions at thermalization which are used to start off the hydrodynamic 
part of the collision fireball evolution. Towards the end of the expansion, 
local thermalization again breaks down because the matter becomes dilute 
and the mean free paths become large. Hence the hydrodynamic evolution 
has to be cut off by hand by implementing a "freeze-out criterion" which 
is discussed at the end of Section [3 

Phenomenological aspects of the expansion are studied in Section [3] 
First we elaborate on central collisions and their characteristics - their 
cooling behavior, their dilution, their expansion rates, etc. We then pro- 
ceed to discuss the special features and possibilities offered by non-central 
collisions, due to the breaking of azimuthal symmetry. We will see how the 
initial spatial anisotropy of the nuclear reaction zone is rapidly degraded 
and transferred to momentum space. This results in a strong momentum 
anisotropy which is easily observable in the measured momentum spectra 
of the finally emitted hadrons. 

Experimental observables reflecting the hydrodynamic fireball history 
are the subject of Section^ These include both momentum and coordinate 
space observables. Momentum space features are discussed in Section ETT1 
We begin by analyzing the angle-averaged transverse momentum spectra of 
a variety of different hadron species for evidence of azimuthally symmetric 
radial flow. The measured centrality dependences of the charged multi- 
plicity near midrapidity, of the mean transverse energy per particle, and 
of the shapes and mean transverse momenta of identified hadron spectra 
are compared with hydrodynamic calculations. We then discuss azimuthal 
momentum anisotropies by decomposing the same spectra into a Fourier 
series with respect to the azimuthal emission angle ip p relative to the reac- 
tion plane. We investigate in particular the mass and centrality dependence 
of the second Fourier coefficient VzIpt), the differential elliptic flow. Com- 
bining the experimental observations with a simple and quite general theo- 
retical arguments, we will build a compelling case for the necessity of rapid 
thermalization and strong rescattering at early times. We will show that 
this provides a very strong argument for the creation of a well-developed 
quark-gluon plasma at RHIC, with a significant lifetime of about 5-7 fm/c 
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Hydrodynamic description of ultrarelativistic heavy-ion collisions 5 

and an initial energy density which exceeds the critical value for color de- 
confinement by at least an order of magnitude. 

In Section f4. 21 we describe how two-particle momentum correlations can 
be used to explore the spatial geometry of the collision fireball at the time of 
hadron emission. This so-called Hanbury Brown- Twiss (HBT) intensity in- 
tcrfcromctry supplements the momentum space information extracted from 
the single particle spectra with information about the size and shape of the 
fireball in coordinate space. Whereas the momentum anisotropics are fixed 
early in the collision, spatial anisotropics continue to evolve until the very 
end of the collision, due to the early established anisotropic flow. By com- 
bining information on the momentum and spatial anisotropics one can hope 
to constrain the total time between nuclear impact and hadronic freeze-out 
rather independently from detailed theoretical arguments. We will show 
that, contrary to the momentum spectra, the experimentally measured 
HBT size parameters are not very well described by the hydrodynamical 
model. Since, contrary to the momentum anisotropics, the HBT correla- 
tions are only fixed at the point of hadronic decoupling, one might suspect 
that they are particularly sensitive to the drastic and somewhat unrealis- 
tic sharp freeze-out criterion empl oyed i n the hydrodynamic simulations. 
However, this so-called HBT-puzzl e^EM [ s shared by most other available 
dynamical models, including those which start hydrodynamically but then 
describe hadron freeze-out kineticallyj ^^lTllS l l^ l S ni and still awaits its 
resolution. 

In our concluding Sectional we give a summary and also highlight the 
need for future studies of uranium- uranium collisions. Since uranium nuclei 
in their ground state are significantly deformed, the long axis being almost 
30% larger than the short one, they offer a significantly deforme d initial 
geometry even for central collisions with complete nuclear overlap 01^221 
The resulting deformed fireballs are much larger and significantly denser 
than those from equivalent semicentral gold-gold collisions, providing bet- 
ter conditions for local thermalization and the validity of hydrodynamic 
concepts even at lower beam energies where data from Au+Au and Pb+Pb 
collisions indicate a gradual breakdown of ideal hydrodynamics!^ 1 Central 
U+U collisions may thus provide a chance to explore the hydrodynamic 
behavior of elliptic flow down to lower collision energies and confirm the 
hydrodynamic predictiorP^^of a non-monotonic structure in its excitation 
function which can be directly related to the quark-hadron phase transition 
and its softening effects on the nuclear equation of state in the transition 
region (see also Section |4~T|| . 
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6 Hydrodynamic description of ultrarelativistic heavy-ion collisions 

2. Formulation of hydrodynamics 
2.1. Hydrodynamic prerequisites 

For macroscopic systems with a large number (say, of the order of Avo- 
gadro's number) of microscopic constituents, thermodynamics takes ad- 
vantage of the fact that fluctuations in the system are small, microscopic 
dynamics drives such systems rapidly to a state of maximum disorder, and 
the system's global behavior can then be expressed in terms of a few macro- 
scopic thermodynamic fields. Thcrmalization happens locally and on mi- 
croscopic time scales which are many orders of magnitude smaller than 
the macroscopic time scales related to the reaction of the system to small 
inhomogeneities of the density, pressure, temperature, etc. Under such con- 
ditions, the system can be described as an ideal fluid which reacts instan- 
taneously to any changes of the local macroscopic fields, by readjusting the 
slope of its particles' momentum distribution (i.e. its temperature) locally 
on an infinitcsimally short time scale. The resulting equations of motion 
for the macroscopic thermodynamic fields are the equations of ideal (i.e. 
non-viscous) hydrodynamics, i.e. the Euler equations and their relativis- 
tic generalizations!^ They describe how macroscopic pressure gradients 
generate collective flow of the matter, subject to the constraints of local 
conservation of energy, momentum, and conserved charges. 

The systems produced in the collision of two large nuclei are much 
smaller: In central Au+Au or Pb+Pb collisions at RHIC energies (i.e. up 
to 200 GeV per nucleon pair in the center of mass system), about 400 nucle- 
ons collide with each other, producing several thousand secondary particles. 
Recent experiments with trapped cold fermionic atoms with tunable inter- 
action strength have shown that systems involving a few hundred thousand 
particles behave hydrodynamically if the local re-equilibration rates are suf- 
ficiently large! I * Similar experiments involving much smaller numbers of 
atoms are under way. In fact, one can argue that the number of particles is 
not an essential parameter for the validity of the hydrodynamic approach, 
and that hydrodynamics does not even rely on the applicability of a particle 
description for the expanding system at all. The only requirement for its 
validity are sufficiently large momentum transfer rates on the microscopic 
level such that relaxation to a local thermal equilibrium configuration hap- 
pens fast on macroscopic time scales. Local thermal equilibrium can also 
be formulated for quantum field theoretical systems which are too hot and 
dense to allow for a particle description because large scattering rates never 
let any of the particles go on-shell. 
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Hydrodynamic description of ultrarelativistic heavy-ion collisions 7 

If the fireballs formed in heavy ion reactions were not expanding, the 
typical macroscopic time scales would be given by the spatial dimensions 
of the reaction zone divided by the speed of sound ~ c/y/3, i.e. of the order 
of lOfm/c. Collective expansion reduces this estimate, and the geometric 
criterion must be replaced by a dynamical one involving the local expansion 
rate ("Hubble constant"), = d^u^ where u^{x) is the local flow 4- 
velocity.™™ 

Typical values for T cxp are of the order of only one 
to several fm/cE^ The hydrodynamic description of heavy- ion collisions 
thus relics on local relaxation times below lfm/c which, until recently, 
was thought to be very difficult to achieve in heavy-ion collisions, causing 
widespread skepticism towards the hydrodynamic approach. The new RHIC 
data have helped to overcome this skepticism, leaving us with the problem 
to theoretically explain the microscopic mechanisms behind the observed 
fast thermalization rates. 

From these considerations it is clear that in heavy-ion collisions a hy- 
drodynamic description can only be valid during a finite interval between 
thermalization and freeze-out. Hydrodynamics can never be expected to 
describe the earliest stage of the collision, just after nuclear impact, dur- 
ing which some of the initial coherent motion along the beam direction 
is redirected into the transverse directions and randomized. The results 
of this process enter the hydrodynamic description through initial condi- 
tions for the starting time of the hydrodynamic stage and for the relevant 
macroscopic density distributions at that time. The hydrodynamic evolu- 
tion is ended by implementing a freeze- out condition which describes the 
breakdown of local equilibrium due to decreasing local thermalization rates. 
These initial and final conditions are crucial components of the hydrody- 
namic model which must be considered carefully if one wants to obtain 
phcnomcnologically relevant results. 

2.2. Hydrodynamic equations of motion 

The energy mom ent um tensor of a thermalized fluid cell in its local rest 
frame is given by-^l T^ st (x) = diag(e(x),p(x),p(x),p(x)) where x labels 
the position of the fluid cell and e(x) and p(x) are its energy density and 
pressure. If in a global reference frame this fluid cell moves with four- velocity 
U^(x) (where = 7 (1, v x , v y , v z ) with 7 = l/yl — v 2 and u^u^ = 1), a 
corresponding boost of T^ st yields the fluid's energy momentum tensor in 
the global frame: 

T^(x)= (e(x) + P {xfju^{x)u u {x) -p{x)g» v . (1) 
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8 Hydrodynamic description of ultrarelativistic heavy-ion collisions 

Note that this form depends on local thermal equilibrium at each point of 
the fluid in its local rest frame, i.e. it corresponds to an ideal fluid where 
dissipative effects can be neglected. The local conservation of energy and 
momentum can be expressed by 

d,T^(x)=0, (z/ = 0, . . . , 3) . (2) 

If the fluid carries conserved charges Ni, with charge densities m(x) in 
the local rest frame and corresponding charge current densities j^{x) — 
ni{x)u^(x) in the global reference frame, local charge conservation is ex- 
pressed by 

d fl j?(x)=0, i = l,...,M. (3) 

Examples for such conserved charges are the net baryon number, electric 
charge, and net strangeness of the collision fireball. 

If local relaxation rates are not fast enough to ensure almost instanta- 
neous local thermalization, the expressions for the energy momentum tensor 
and charge current densities must be generalized by including dissipative 
terms proportional to the tran sport c oefficients for diffusion, heat conduc- 
tion, bulk and shear viscosity! ' The solution of the correspondingly 
modified equations is very challenging!^ We will later discuss some first 
order viscous corrections in connection with experimental observables. 

2.3. The nuclear equation of state 

The set ll'-'l.'H) of 4 + M differential equations involves 5 + M undetermined 
fields: the 3 independent components of the flow velocity, the energy density, 
the pressure, and the M conserved charge densities. To close this set of 
equations we must provide a nuclear equation of state p(e, rii) which relates 
the local thermodynamic quantities. We consider only systems with zero 
net strangeness and do not take into account any const rain ts from charge 
conservation which are known to have only minor effects!^ This leaves the 
net baryon number density n as the only conserved charge density to be 
evolved dynamically. 

The equation of state for dense systems of strongly interacting particles 
can either be modeled or extracted from lattice QCD calculations. We use 
a combination of t hes e two possibilities: In the low temperature regime, 
we follow HagedorrP^ and describe nuclear matter as a noninteracting gas 
of hadronic resonances, summing over all experimentally identified"^ res- 
onance states ElsiSI As the temperature is increased, a larger and larger 
fraction of the available energy goes into the excitation of more and heavier 
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resonances. This results in a relatively soft equation of state ("EOS H") 
with a smallish speed of sound: c 2 = dp/de ~ 0.15P3I 

As the available volume is filled up with resonances, the system ap- 
proaches a phase transition in which the hadrons overlap and the micro- 
scopic degrees of freedom change from hadrons to deconfined quarks and 
gluons. Due to the large number of internal quark and gluon degrees of 
freedom (color, spin, and flavor) and their small or vanishing masses, this 
transition is accompanied by a rather sudden increase of the entropy density 
at a critical temperature T cr i t . Above the transition, the system is modeled 
as a noninteracting gas of massless u, d, s quarks and gluons, subject to an 
external bag pressure jB.ESlThe corresponding equation of state p = |e— |i? 
is quite stiff and yields a squared sound velocity c 2 — dp/de — -| which is 
more than twice that of the hadron resonance gas. In the following we refer 
to this equation of state as "EOS I" . 




e (GeV/fm J ) 



Fig. 1. Equation of state of the Hagedorn resonance gas (EOS H), an ideal gas of mass- 
less particles (EOS I) and the Maxwellian connection of those two as discussed in the 
text (EOS Q). The figure shows the pressure as function or energy density at vanishing 
net baryon density. 



We match the two equations of state by a Maxwell construction, ad- 
justing the bag constant B X I A = 230 MeV such that for a system with zero 
net bar yon density the transition temperature coincides with lattice QCD 
results PEI We choosd^l T cr it = 164 MeV and call the resulting combined 
equation of state "EOS Q" . It is plotted as p(e) at vanishing net baryon 
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10 Hydrodynamic description of ultrarelativistic heavy-ion collisions 

density and strangeness in Figure together with the hadron resonance 
gas EOS H and the ideal gas of massless partons, EOS I. The Maxwell con- 
struction inevitably leads to a strong first order transition, with a large 
latent heat Aei a t = 1.15 GeV/fm 3 (between upper and lower critical val- 
ues for the energy density of eq = 1.6 GeV/fm 3 and en = 0.45 GeV/fm 3 )!^ 
This contradicts lattice QCD resultJ^El which suggest a smoother transi- 
tion (either very weakly first order or a smooth cross-over). However, the 
total incre ase o f the entropy density across the transition as observed in the 
lattice datcPEl i s well reproduced by the model, and it is unlikely that the 
artificial sharpening of the transition by the Maxwell construction leads to 
significant dynamical effects. This is in particular true since the numerical 
algorithm used to solve the hydrodynamic equations tends to soften shock 
discontinuities such as those which might be generated by a strong first 
order phase transition. We will return to the possible influence of details of 
the EOS on certain observables in Section 0] when we discuss experimental 
data. 

2.4. Initialization 

As discussed in Section 12.11 the initial thermalization stage in a heavy- 
ion collision lies outside the domain of applicability of the hydrodynamic 
approach and must be replaced by initial conditions for the hydrodynamic 
evolution. Different authors have explored a variety of routes to arrive at 
reasonable such initial conditions. For example, one can try to treat the two 
colliding nuclei as two interpenetrating cold fluids feeding a third hot fluid 
in the reaction center ("three-fluid dynamics'^} . This requires modelling 
the source and loss terms describing the exchange of energy, momentum 
and baryon number among the fluids. Alternatively, one can model the 
early stage kinetically, using a transport model such as the parton cascades 
VNPHI VNI/BMS22 MPCfEl AMPT^or one of several other available 
transport codes to est imate the initial energy and entropy distributions in 
the collision regional before switching to a hydrodynamic evolution. 

However, the microscopic effects which generate the initial entropy are 
still poorly understood, and it is quite likely that, due to the high density 
and collision rates, transport codes which simulate the solution of a Boltz- 
mann equation using on-shell particles are not really valid during the early 
thermalization stage. In our own calculations, we therefore simply parame- 
terize the initial transverse entropy or energy density profile geometrically, 
using an optical Glauber model calculation^ to estimate the density of 
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Hydrodynamic description of ultrarelativistic heavy-ion collisions 11 

wounded nucleons and binary nucleon-nucleon or parton-parton collisions 
in the plane transverse to the beam and superimposing a "soft" component 
(scaling with the number of wounded nucleons) and a "ha rd" com ponent 
(scaling with the number of binary collisions) in such that the 

experiment ally ob served rapidity density of charged hadrons at the end of 
the collisiorPHEi and its dependence on the collision centralit)!^^ are 
reproduced.^ 

For the Glauber calculation we describe the density distributions of the 
colliding nuclei (with mass numbers A and B) by Woods-Saxon profiles, 

P^) = e (r-«ff/€ + l ' W 

with the nuclear radius R A = {1.12 A 1 / 3 — 0.86 A -1 / 3 ) fm and the surface 
diffuscness £ = 0.54fmp21 The nuclear thickness function is given by the 
optical path-length through the nucleus along the beam direction: 

/oo 
dzp A (x,y, z). (5) 
-oo 

The coordinates x, y parametrize the transverse plane, with x pointing in 
the direction of the impact parameter b (such that (x, z) span the reaction 
plane) and y perpendicular to the reaction plane. For a non-central collision 
with impact parameter b, the density of binary nucleon-nucleon collisions 
n BC at a point (x,y) in the transverse plane is proportional to the product 
of the two nuclear thickness functions, transversally displaced by b: 

n BC (x,y;b) = <r T A (x + b/2, y) T B (x - 6/2, y). (6) 

(To is the total inelastic nucleon-nucleon cross section; it enters here only as 
a multiplicative factor which is later absorbed in the proportionality con- 
stant bet ween n BC (x, y; b) and the "hard" component of the initial entropy 
deposition!^ Integration over the transverse plane (the (x, y)-plane) yields 
the total number of binary collisions, 

N vc( b )= / dxdy n BC {x,y; b) (7) 



Its impact parameter dependence, as well as that of the maximum density 
of binary collisions in the center of the reaction zone, n BC (0, 0; 6), are shown 
as the dashed lines in Fig. 

The "soft" part of the initial entropy deposition is assumed to scale 
with the density of "wounded nucleons" j™ defined as those nucleons in the 
projectile and target which participate in the particle production process 
by suffering at least one collision with a nucleon from the other nucleus. 
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The Glauber model gives the following transverse density distribution of 
wounded nucleons!^ 



n wn {x,y\b) = T A (x + b/2,y) 1 - ( 1 - 



T B (x-b/2,y) 1-1- 



<jqTb(x - b/2, yY 
B 

<7gT A (x + b/2,y )' 
A 



(8) 



Here the value Co of the total inelastic nucleon-nucleon cross section plays 
a more important role since it influences the sh ape of the transverse den- 
sity distribution n WN (a;, y; b), and its dependence^ on the collision energy 
y/s must be taken into account. The total number of wounded nucleons is 
obtained by integrating Eq. JSJ over the transverse plane. Its impact pa- 
rameter dependence, as well as that of the maximum density of wounded 
nucleons in the center of the reaction zone, n WN (0,0;6), are shown as the 
solid lines in Fig. 
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Fig. 2. Left: Number of wounded nucleons and binary collisions as a function of impact 
parameter, for Au+Au collisions y/s = 130AGeV (cro = 40 mb). Right: Density of 
wounded nucleons and binary collisions in the center of the collision as a function of 
impact parameter. 



Our hydrodynamic calculations were done with initial conditions which 
ascribed 75% of the initial entropy production to "soft" processes scaling 
with n WN (x, y; b) a nd 25% to "hard" processes scaling with n BC (x .y; b) . 
This was found^l to give a reasonable description of the measure d ' " 
centrality dependence of the produced charged particle rapidity density per 
participating ("wounded") nucleon. 

Figure [3 shows that exploring the centrality dependence of heavy-ion 
collisions provides access to rich physics: With increasing impact parameter 
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Hydrodynamic description of ultrarelativistic heavy-ion collisions 13 

both the number of participating nucleons and the volume of the created 
fireball decrease. Except for effects related to the deformation of the reac- 
tion zone in non-central collisions, increasing the impact parameter is thus 
equivalent to colliding smaller nuclei, eventually reaching the limit of pp 
collisions in the most peripheral nuclear collisions. Furthermore, at fixed 
beam energy, the initially deposited entropy and energy densities decrease 
with increasing impact parameter. To a limited extent, heavy-ion collisions 
at fixed beam energy but varying impact parameter are therefore equiv- 
alent to central heavy-ion collisions at varying beam energy, i.e. one can 
map sections of the "excitation function" of physical observables without 
changing the collision energy, but only the collision centrality. 

In one respect, however, non-central collisions of large nuclei such as 
Au+Au are fundamentally different from central collisions between lighter 
nuclei: A finite impact parameter breaks the azimuthal symmetry inherent 
in all central collisions. In a strongly interacting fireball, the initial geomet- 
ric anisotropy of the reaction zone gets transferred onto the final momentum 
spectra and thus becomes experimentally accessible. As we will see, this 
provides a window into the very early collision stages which is completely 
closed in central collisions between spherical nuclei. (The same information 
is, however, accessible, with even better statistics due to the larger overlap 
volume and number of produced particles, in completely central collisions 
between deformed nuclei, such as W+W or U+U.) 



10 




-10 -5 5 10 2 4 6 8 10 12 



x (fm) b (fm) 

Fig. 3. Density of binary collisions in the transverse plane for a Au+Au collision with 
impact parameter 6 = 7 fm (left). Shown are contours of constant density together with 
the projection of the initial nuclei (dashed lines). The right plot shows the geometric 
eccentricity as a function of the impact parameter for the wounded nucleon and binary 
collision distributions. 
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14 Hydrodynamic description of ultrarelativistic heavy-ion collisions 

The left panel of Fig. [21 shows the distribution of binary collisions in 
the transverse plane for Au+Au collisions at impact parameter 6=7 fm. 
Shown are lines of constant density at 5, 15, 25, ... % of the maximum 
value. The dashed lines indicate the Woods-Saxon circumferences of the 
two colliding nuclei, displaced by ±6/2 from the origin. The obvious geo- 
metric deformation of the overlap region can be quantified by the spatial 
eccentricity 

(v 2 - t 2 ) 

where the averages are taken with respect to the underlying density (n WN 
or n BC or a combination thereof, depending on the exact parametrization 
used). The centrality dependence of e x is displayed in the right panel of 
Fig. El 

2.5. Decoupling and freeze-out 

As already mentioned in Section 12. II the hydrodynamic description begins 
to break down again once the transverse expansion becomes so rapid and 
the matter density so dilute that lo cal thermal equilibrium can no longer be 
maintained. Detailed studieS^LHI comparing local mean free paths with the 
overall size of the expanding fireball and the local Hubble radius (inverse 
expansion rate) have shown that bulk freeze-out happens dynamically, i.e. 
it is driven by the expansion of the fireball and not primarily by its finite 
size. This is similar to the decoupling of the primordial nuclear abundances 
and the cosmic microwave background in the early universe which was also 
entirely controlled by the cosmic expansion rate. Nonetheless, some part of 
the initially produced matter never becomes part of the hydrodynamic fluid, 
but decouples right away even though no transverse flow has developed yet. 
Figure shows that already at initialization the density distribution has 
dilute tails where the mean free path is never short enough to justify a 
hydrodynamic treatment. These tails should be considered as immediately 
frozen out, i.e. they describe particles which carry their momentum infor- 
mation directly and without further modification to the detector. Their 
decoupling is obviously not a result of (transverse) dynamics, but a geo- 
metric effect. However, for both geometric and dynamical freeze-out the 
local scattering rate (density times cross section) is the controlling factor, 
with the density showing the largest variations across the fireball, and it 
was founcP^ESI ^ na |. ^ e hypersurface along which the local mean free path 
begins to exceed the local Hubble radius or the global fireball size can be 
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characterized in good approximation as a surface of constant temperature. 
Note that for almost baryon-antibaryon symmetric systems such as the ones 
generated near midrapidity at RHIC, the entropy density, energy density, 
particle density and temperature profiles are directly related and all have 
similar shapes. A surface of constant temperature is therefore, in excellent 
approximation, also a surface of constant energy and particle density. 

A traditional way of describing the breakdown of hydrodynamics and 
particle freeze-out is the Cooper-Frye prescriptiorl^ which postulates a 
sudden transition from perfect local thermal equilibrium to free streaming 
of all particles in a given fluid cell once the kinetic freeze-out criterion (ob- 
tained, for example, in the way just described, by comparing local scattering 
and expansion rates etc.) is satisfied in that cell. In the Cooper-Frye formal- 
ism, one first lets hydrodynamics run up to large times, then determines the 
space-time hypersurface along which the hydrodynamic fluid cells first 
pass the freeze-out criterion, and computes the final spectrum of particles 
of type i from the formula 

e W- = h 7\i = WT* / P (10) 

d 6 p dypTdpTdipp (2-7r) d J E 

where drajix) is the outward normal vector on the freeze-out surface £(x) 
such that p^fi d 3 <Tfj, is the local flux of particles with momentum p through 
this surface. For the phase-space distribution / in this formula one takes 
the local equilibrium distribution just before decoupling, 

boosted with the local flow velocity u M (x) to the global reference frame by 
the substitution E — > p ■ u(x). /itj(x) and T(x) are the chemical potential 
of particle species i and the local temperature along E, respectively. The 
temperature and chemical potentials on £ are computed from the hydrody- 
namic output, i.e. the energy density e, net baryon density n and pressure 
p, with the help of the equation of stateP^ 

This formalism is used to calculate the momentum distribution of all di- 
rectly emitted hadrons of all masses. Unstable resonances are then allowed 
to decay (we include all strong decays, but consider weakly decaying par- 
ticles as stable), taki ng th e appropriate branching ratio of different decay- 
channels into account!"^ The contribution of the decay products is added 
to the thermal momentum spectra of the directly emitted stable hadrons 
to give the total measured particle spectra!^ 
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Initial particle production at high px is controlled by hard QCD pro- 
cesses which produce transverse momentum spectra which strongly deviate 
from an exponential thermal or hydrodynamic shape. Since the momen- 
tum transfer per collision is limited, such particles require a larger number 
of rescatterings than soft particles for reaching thermal equilibrium. They 
have a much higher chance of escaping from the fireball before being ther- 
malized than soft hadrons. This is not captured by the Cooper-Frye formula 
which freezes out all particles in a given fluid cell together, irrespective of 
their momenta. A modified freeze-out criterium which takes the momen- 
tum dependence of the escape probability into account has recently been 
advocated.EZEHEi We will discuss phenomena at large transverse momenta 
in the last part of Section Em 

The Cooper-Frye formula has another shortcoming which materializes 
if the freeze-out normal vector d 3 a fJ, (x) is spacelike (as it happens in cer- 
tain regions of our hydrodynamic freeze-out surfaces), in which case the 
Cooper-Frye integral also counts (with negative flux) particles entering the 
thermalized fluid from outside. However, simple attempts to cut these con- 
tributions by hand^l generate problems with energy-momentum conser- 
vation, and only recently a possible resolution of this problem has been 
found.EU 

Clearly, any Cooper-Frye like prescription implementing a sudden tran- 
sition from local equilibrium (infinite scattering cross section) to free- 
streaming (zero cross section) is ultimately unrealistic and should be re- 
placed by a more realistic prescription. A preferred procedure would be 
the transition from hydrodynamics t o kinetic transport theory just before 
hydrodynamics begins to break down,^ ^^^^^ thus allowing for a 
gradual decoupling process which is fully consistent with the underlying 
microscopic physics. Such a realistic treatment of the freeze-out process is 
clearly much more involved than the Cooper-Frye formalism, and so far it 
has not led to strong qualitative differences for the emitted hadron spec- 
tra, even though in detail some phenomenological advantages of the hybrid 
(hydro+transport) approach can be identified Also, the crucial question 
whether for rapidly expanding heavy-ion fireballs there is really an over- 
lap window where both the macroscopic hydrodynamic and the microsopic 
transport approach using on-shell particles work simultaneously has not 
been finally settled. Most of the results presented in this review have there- 
fore been obtained using the simple Cooper-Frye freeze-out algorithm. 
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2.6. Longitudinal boost invariance 

Most of the observables to be discussed below have been collected near 
mid-rapidity. This region is of particular interest as one expects there the 
energy- and particle density to be the largest, giving the clearest signals of 
the anticipated phase-transition, and many components of the large heavy- 
ion experiments have therefore been optimized for the detection of midra- 
pidity particles. Furthermore, rapidity distributions are more difficult to 
analyze theoretically than transverse momentum distributions since they 
are strongly affected by "memories of the pre-collision state" : Whereas all 
transverse momenta are generated by the collision itself, a largely unknown 
fraction of the beam-component of the momenta of the produced hadrons 
is due to the initial longitudinal motion of the colliding nuclei. In hydro- 
dynamics one finds that final rapidity distributions are very sensitive to 
the initialization along the beam direction, and that hydrod yna mic evolu- 
tion is not very efficient in changing the initial distributions!^ Collective 
transverse effects are thus a cleaner signature of the reaction dynamics than 
longitudinal momentum distributions, and the best way to isolate oneself 
longitudinally from remnants of the initially colliding nuclei is by going 
as far away as possible from the projectile and target rapidities, i.e. by 
studying midrapidity hadron production. 

Near midrapidity one is far from the kinematic limits imposed by the 
finite collision energy, and the microscopic processes responsible for parti- 
cle production, scattering, thermalization and expansion should therefore 
be locally the sam e ev erywhere and invariant under limited boosts along 
the beam direction!^ In a hydrodynamic description this implies a boost- 
invariant longitudinal flow velocitj^l whose form is independent of the 
transverse expansion of the fluid while the latter is identical for all trans- 
verse planes in their respective longitudinal rest frames. Under these as- 
sumptions the analytically solvable longitudinal evolution decouples from 
the transverse evolutio n"^ ' u l which greatly reduces the numerical task of 
solving the hydrodynamic equations of motion. Limitations of Bjorken's 
solution and boost invariance will be discussed in Section T4. 1.41 but most 
of the review reports results obtained under the assumption of longitudinal 
boost invariance. 



Bjorkcn showed^ that the boost invariant longitudinal flow field has 
the scaling (Hubble) form v z — z/t and that the hydrodynamic equa- 
tions preserve this form in proper time if the initial conditions for the 
thermodynamic variables do not depend on space-time rapidity rj — 
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\ ln[(t+z) / (t—z)]. With this profile the flow 4- velocity can be parametrized 
as u M = cosh y± (cosh 77, v x , v yi sinh^), where the transverse flow rapidity 
yj_(r,r) does not depend on 77 and is related to the radial flow velocity 
v r = ^J v x^ v y °^ midrapidity rj = by v r (r, r, 77=0) = tanh yj_(r, r). It is 
then sufficient to solve the hydrodynamic equations for v r at z = 0, and 
the transverse velocity at other longitudinal positions is given bjEIl 

v r (t,r,z) = v r (T,r,r 1 ) = Vr{T >l> 0) . (12) 

cosn?7 

Solving the hydrodynamic equations in the transverse plane r = (x, y) 
with longitudinally boost-invariant boundary conditions becomes easiest 
after a coordinate transformation from (z,t) to longitudinal proper time r 
and space-time rapidity 77: 



(t, x, y, z) — > x m = (r, x, y, 77) 



t = t cosh 77 r = vt 2 — z 2 

(13) 

z = Tsinh?7 77 = Artanh(z/t) . 

In these coordinates, the equations of motions becom^ 
T T \ T + {vxT TT ) x + {v y T TT ) y = -i (T TT + p) - {pv x ), x - (pv y ), v , (14) 
T TX , T + (v x T TX ) x + (v v T™) = - p>x -It™, (15) 

. y T 

T T \ T + (y x T Ty ) x + (v y T T y) = -p y -- T Ty , (16) 

' y T 

^P-n = 0, (17) 

f,r + (V x f), x + {Vyf),y = -- f , (18) 

T 

where the lower case comma indicates a partial derivative with respect to 
the coordinate following it. One sees that the evolution in 77-direction is 
now trivial, and that only 4 coupled equations remain to be solved. 
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3. Phenomenology of the transverse expansion 

In this section we study the transverse fireball expansion at midrapidity as 
it follows from the hydrodynamic equations of motion fSections l2.2l and l2.ti|) 
with the equation of state described in Section f2.3l and the initial conditions 
from Section EH 

In the first part of this section we study central collisions (6 = 0). These 
are used to tune the initial conditions of the calculations, by requiring the 
calculation to reproduce the measured rapidity density of charged hadrons 
at midrapidity and the shape of the pion and proton spectra in central 
collisions. For Au+Au collisions with a center of mass energy of 200 GeV per 
nucleon pair, we find for the initial equilibration time (i.e. for the beginning 
of the hydrodynamic stage) r oqu = 0.6 fm/c and an initial entropy density in 
the center of the fireball of s cqu = 110 fm -3 !^ Freeze-out occurs when the 
energy density drops below edcc = 0.075 GeV/fm 3 . How these parameters 
are fixed will be described in some detail in Section 14.11 In the present 
section we will simply use them to illustrate some of the characteristic 
features of the transverse hydrodynamic expansion. 

In the second part we address non-central collisions and discuss the spe- 
cial opportunities provided by the breaking of azimuthal symmetry in this 
case. We discuss how the initial spatial deformation transforms rapidly into 
a momentum space anisotropy which ultimately manifests itself through a 
dependence of the emitted hadron spectra and their momentum correlations 
on the azimuthal emission angle relative to the reaction plane. 

3.1. Radial expansion in central collisions 

As seen from the terms on the right hand side in Eqs. 114|I - (|16J1 . the driv- 
ing force for the hydrodynamic expansion are the transverse pressure gra- 
dients which accelerate the fireball matter radially outward, building up 
collective transverse flow. As a result, the initial one-dimensional boost- 
invariant expansion along the beam direction gradually becomes fully three- 
dimensional. For the adiabatic (ideal fluid) expansion discussed here, this 
implies that the entropy and other conserved charges spread out over a vol- 
ume which initially grows linearly with time, but as time evolves increases 
faster and ultimately as r 3 . Accordingly, the entropy and baryon densi- 
ties follow inverse power laws ~ r~ a with a "local expansion coefficient" 
a = — | | t " - which changes smoothly from 1 to 3. 

This is seen in the left panel of Figure 0] which shows a double logarith- 
mic plot of the entropy density at three points in the fireball (at the origin 
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Fig. 4. Left panel: Time evolution of the entropy density at three different points in 
the fireball (0, 3, and 5 fm from the center). Dashed lines indicate the expectations for 
pure one-dimensional and three-dimensional dilution, respectively. Right panel: Time 
evolution of the temperature at the same points. The plateau at T = f64MeV results 
from the transition of the corresponding fluid cells through the mixed phase. 

R = as well as 3 and 5 fm away from it) as a function of time. At early 
and late times the solid lines representing the numerical solution are seen 
to follow simple r _1 and r~ 3 scaling laws, indicated by dashed lines. For 
an ideal gas of massless particles (such as our model for the quark-gluon 
plasma above the hadronization phase transition) s ~ T 3 , and ar" 1 scaling 
of the entropy density translates into a r -1 / 3 scaling of the temperature. 
The right panel of Figure 0] shows that the numerical solution follows this 
simple scaling rather accurately almost down to the phase transition tem- 
perature T cr i t . At this point the selected fluid cell enters the mixed phase 
and the temperature remains constant until the continued expansion di- 
lutes the energy density below the lower critical value en = 0.45 GeV/fm 3 . 
This is an artefact of the Maxwell construction employed in Section E"3l for 
a more realistic rapid but smooth crossover between the QGP and hadron 
gas phases the horizontal plateau in Figure0]would be replaced by a similar 
one with smoothed edges and a small but finite slope. As the fluid cell exits 
the phase transition on the hadronic side, its temperature is seen to drop 
very steeply; this is caused not only by the now much more rapid three- 
dimensional expansion, but also by the different temperature dependence of 
s(T) in the hadronic phase, generated by the exponential dependence of the 
phase-space occupancy on the hadron rest masses. The horizontal dashed 
line indicates the freeze-out temperature of about 100 MeV (see Section 
14.1(1 : one sees that in central collisions freeze-out occurs about 15 fm/c 
after equilibration. 

Figure compares the time-dependence of the local expansion coeffi- 
cient a = — = — 7 §^ with that of the local expansion rate divided by 



5, 2008 7:22 



WSPC/Trim Size: 9in x 6in for Review Volume 



Hydrodynamic description of ultrarelativistic heavy-ion collisions 



21 



the boost-invariant longitudinal expansion rate, T(d fl u fJ, )W^The horizontal 
dashed lines indicate the expectations for pure 1-dimensional longitudinal 
(Bjorken-like) expansion (a = r(3^ii A1 ) = 1) and for 3-dimensional isotropic 
radial (Hubble-like) expansion (a = r(d^u^) = 3), respectively. One sees 
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Fig. 5. Left panel: Time evolution of the local expansion coefficient a = — 9(ln s) /<9(ln r) 
at three different fireball locations. The horizontal dashed lines indicate expectations for 
pure Bjorken (a = 1) and Hubble-like expansions (a = 3). Right panel: The local 
expansion rate 9 m m m multiplied by time, again compared to Bjorken and Hubble-like 
scaling expansions. Note that the expansion coefficient a decreases with radial distance 
from the center whereas the expansion rate shows the opposite behavior. 



that in the numerical solutions both quantities increase with time from 
a Bjorken-like behavior initially to a Hubble-like behavior at later times. 
Weak structures in the time evolution at the beginning and end of the 
mixed phase are probably due to our artificially sharp phase transition and 
should disappear for a realistic equation of state. Note that both the lo- 
cal expansion coefficient a and the normalized local expansion rate r d-u 
exceed the limiting global value of 3 at large times. This does not violate 
causality, bu t is due to the existence of density gradients and their time 
evolution!^ 

At late times the expansion rate r d-u is larger for points at the edge of 
the fireball than in the center, again due to the stronger density and pres- 
sure gradients near the edge. In contrast, the local dilution rate a shows 
the opposite dependence on the radial distance, being smaller at large ra- 
dial distance than in the center. This reflects the transport of mat ter from 
the center to the edge, due to radial flow and density gradients!^ The 
relation between a and r d-u can be established by using entropy con- 
servation, 9 m (sm m ) = 0, to write d-u= — (u-ds)/s. Assuming longitudinal 
boost-invariance and a temporal power law s(r, r) = so(r)(^ L ) Q for the lo- 
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cal entropy dilution rate one finds the relational 

a„«* = 7 ("-«r^). (19) 

The last term involving the radial flow and radial entropy density gradient 
is positive, especially at large radii, explaining the different ordering of the 
three curves in the left and right panels of Figure El 

To further illustrate the transverse dynamics we show in Figure El the 
radial velocity function of time and radial distance. The left panel 

shows the time evolution of v r at fixed radii of 3 and 5 fm. After a steep 
initial rise of v r , the radial velocity at fixed position R is seen to decrease 
with time while the matter there changes from QGP into a hadron gas. 
Inside the mixed phase the pressure is constant (i.e. the pressure gradient 
vanishes) and the matter is not further accelerated. As a result, the system 
expands without acceleration, with rapidly flowing matter moving to larger 
radii while more slowly moving matter from the interior arrives at the 
fixed radius R. Only after the mixed phase has completely passed through 
the radius R does the radial expansion accelerate again, caused by the 
reappearance of pressure gradients. 




048 12 16 02468 
x (fm/c) r (fm) 



Fig. 6. Left panel: Time evolution of the collective radial velocity at 3 and 5 fm radial 
distance from the origin. The short vertical lines in the figure indicate the beginning 
and end of the hadronization phase transition and the freeze-out time. Right panel: The 
radial velocity profile v r (r) at four different times. The dashed line indicates a linear 
profile v r (r) = £r with slope £ = 0.07 fm -1 . 

The right panel of Figure shows the radial velocity profile at selected 
times. Initially v r = whereas the pressure profile p(r) features strong radial 
gradients (especially near the surface), except for a moderately thin layer 
between the QGP and hadron gas phases where the matter is in the phase 
transition and the corresponding softness of the equation of state does not 
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allow for pressure gradients. Accordingly, the initial acceleration happens 
mostly in the outer part of the QGP core and in the hadron gas shell, 
with no acceleration in the mixed phase layer which only gets squashed by 
the accelerating matter pushing out from the interior. This is clearly seen 
in the solid line in Figure which represents the radial velocity profile at 
t = 1 fm/c. As time proceeds the structure caused by the weak acceleration 
in the mixed phase gets washed out and the velocity profile becomes more 
uniform. It very rapidly approaches a nearly linear shape v r (r)^^r with 
an almost time-independent limiting slope of £ ~ 0.07 fm _1 . 




R T (fm) R^fm) 



Fig. 7. The transverse flow rapidity y± as a function of radial distance r along a surface 
of constant energy density e = 0.45 GeV/fm 3 , for Pb+Pb collisions at the SPS (left) and 
for Au+Au collisions at R HIC (right) lial Three different equations of state have been 
explored in this figureJEl with LH8 corresponding most closely to EOS Q shown in 
Figure The dashed and solid line segments subdivide the surface into 5 equal pieces 
of 20% each with respect to the entropy flowing through the surface. 

The near constancy of this slope implies that one also obtains an almost 
linear transverse velocity profile along a hypersurface of fixed temperature 
or energy density rather than fixed time. Figure \7\ compares the radial flow 
rapidity profile y±(r) for Pb+Pb or Au+Au collisions at SPS and RHIC en- 
ergies for three different equations of state, ^ with LH8 corresponding most 
closely to EOS Q shown in Figure ^ Figure provides welcome su pport 
for the phenomenologically very successful blast-wave parametrizatio n" ' 
which is usually employed with a linear transverse velocity or rapidity pro- 
file for reasons of simplicity. (Note that for the range of velocities covered 
in the figure the difference between rapidity y± and velocity v r = tanhy^ 
can be neglected.) 



February 5, 2008 7:22 



WSPC/Trim Size: 9in x 6in for Review Volume 



24 Hydrodynamic description of ultrarelativistic heavy-ion collisions 

As discussed in Sections l2.1l and l2"51 the freeze-out of particle species i is 
mostly controlled by the competition between the macroscopic expansion 
time scal d^ l & LI Texp = (g.-u)- 1 and the microscopic scattering time scale 
r scatt = V J2j ( a ij v ij) Pjl here the sum goes over all particle species j (with 
density pj) in the fireball and (<JijVij) are the corresponding thermally av- 
eraged and velocity weighted s cattering cross sections. The scattering rates 
drop steeply with temperature P^I^I^ I enabling us to idealize freeze-out as 
a relatively sudden process which happens along a freeze-out surface £ of 
approximately constant decoupling temperature Td cc - The most important 
contributions to the local scattering rate arise from 7r+meson (due to their 
large abundance) and 7r+(anti)baryon collisions (due to their large resonant 
cross sections). At RHIC the total baryon density is somewhat lower than 
at the SPS, due to the smaller baryon chemical potential (note that this 
does not reduce the contribution from baryon-antibaryon pairs!), and one 
thus expects that, at the same temperature, the mean scattering time r scat t 
should be slightly longer at RHIC than at the SPS. The magnitude of this 
effect should be small, however, and its sign could even be reversed if at 
RHIC the pion phase-space is significantly oversaturatedP^ 

On the other hand, the expansion time scale r Gxp = (<9-u) _1 does change 
significantly between SPS and RHIC: For boost-invariant longitudinal flow 
and a linear transverse flow rapidity profile y± =£r (as suggested by Figs. El 
and[7|l the expansion rate is easily calculated aa^ 

cosh(£r) / , . sinh(£r)\ 1 ^ . „. 

d-u= K ±-L + £ cosher) + ^ 1 « - + 2£, (20) 

r \ s r / T 

where the approximate expression^! holds in the region (r<l. Equa- 
tion 112011 gives r(d-u) = 1 + 2£r; reading off £ w 0.07 fm _1 from Figures 
and at RHIC energies, this linear function reproduces well the almost lin- 
ear behavior seen in the right panel of Fig. On the other hand, Figured 
shows at SPS energies a transverse flow rapidity slope t hat is only about 2/3 
of the value at RHIC. At freeze-out (r doc ~ 15 - 17fm/P^) the expansion 
rate at RHIC is thus about 25% larger than at the SPS ((<3-it)dec ~ 0.21 fm _1 
for Au+Au at y/s= 130 A GeV vs. {d-u) dcc « O.lBfm- 1 for Pb+Pb at 
yfs = 1 7 A GeV) . The corresponding "Hubble times" at freeze-out are 
r^fp (RHIC) w 4.8 fm/c and T c d x °=(SPS) w 6.1 fm/c. Combining this with the 
already mentioned rather similar microscopic scattering time scales at both 
energies one is led to the conclusion that at RHIC freeze-out should happen 
at a somewhat higher decoupling temperature than at the SPS. 
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3.2. Anisotropic flow in non-central collisions 

In Section TIM we have already addressed some of the great opportunities 
offered by non-central collisions. The most important ones are related to the 
broken azimuthal symmetry, introduced through the spatial deformation of 
the nuclear overlap zone at non-zero impact parameter (see Figure[2J| • If the 
system evolves hydrodynamically, driven by its internal pressure gradients, 
it will expand more strongly in its short direction (i.e. into the direction 
of the impact parameter) than perp endicular to the reaction plane where 
the pressure gradient is smaller.L^ This is shown in Figure [3] where con- 
tours of constant energy density are plotted at times 2, 4, 6 and 8 fm/c 
after thermalization. The figure illustrates qualitatively that, as the system 
evolves, it becomes less and less deformed. In addition, some interesting 
fine structure develops at later times: After about 6 fm/c the energy den- 
sity distribution along the x-axis becomes non-monotonous, forming two 
fragments of a shell that enclose a little 'nut' in the center V-** Unfortun- 
mately, when plotting a cross section of the profiles shown in Figure one 
realizes that this effect is rather subtle, and it was also found to be fragile, 
showing a strong sensitivity to details of the initial density profiled 
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Fig. 8. Contours of constant energy density in the transverse plane at different times 
(2, 4, 6 and 8 fm/c after equil ibrat ion) for a Au+Au collision at ^/snn = 130 GcV 
and impact parameter 6 = 7 fm p l ' ^ ' Contours indicate 5, 15, . . . , 95 % of the maximum 
energy density. Additionally, the black solid, dashed and dashed-dotted lines indicate the 
transition to the mixed-phase, to the resonance gas phase and to the decoupled stage, 
where applicable. 

A more quantitative characterization of the contour plots in Figure [S] 
and their evolution with time is provided by defining the spatial eccentricity 



(y 2 ~x 2 ) 



(21) 



(y 2 + x 2 ) ' 

where the brackets indicate an average over the transverse plane with 
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the local energy density e(x, y;r) as weight function, and the momentum 
anisotropy 

J dxdy (T xx - T yy ) 



e p (r) 



(22) 



/ dxdy (T xx + TVV) ' 

Note that with these sign conventions, the spatial eccentricity is positive 
for out-of-plane elongation (as is the case initially) whereas the momentum 
anisotropy is positive if the preferred flow direction is into the reaction 
plane. 




x (fm/c) 



Fig. 9. Time evolution of the spatial eccentricity e x and the momentum anisotropy e p 
for Au+Au collisions at RHIC with b = 7 fmE3 

Figure shows the time evolution of the spatial and momentum 
anisotropies for Au+Au collisions at impact parameter £> = 7fm, for RHIC 
initial conditions with a realistic equation of state (EOS Q, solid lines) and 
for a much higher initial energy density (initial temperature at the fire- 
ball center = 2 GeV) with a massless ideal gas equation of state (EOS I, 
dashed lines) P*^ The initial spatial asymmetry at this impact parameter is 
ex(Tequ) = 0.27, and obviously e p (T oqu ) = since the fluid is initially at rest 
in the transverse plane. The spatial eccentricity is seen to disappear before 
the fireball matter freezes out, in particular for the case with the very high 
initial temperature (dashed lines) where the source is seen to switch orien- 
tation after about 6 fm/c and becomes in-plane-elongated at late times!^ 
One also sees that the momentum anisotropy e p saturates at about the 
same time when the spatial eccentricity e x vanishes. All of the momentum 
anisotropy is built up during the first 6 fm/c. 
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Near a phase transition (in particular a first order transition) the equa- 
tion of state becomes very soft, and this inhibits the generation of transverse 
flow. This also affects the generation of transverse flow anisotropics as seen 
from the solid curves in Figure The rapid initial rise of e p suddenly stops 
as a significant fraction of the fireball matter enters the mixed phase. It 
then even decreases somewhat as the system expands radially without fur- 
ther acceleration, thereby becoming more isotropic in both coordinate and 
momentum space. Only after the phase transition is complete and pressure 
gradients reappear, the system reacts to the remaining spatial eccentric- 
ity by a slight further increase of the momentum anisotropy. The softness 
of the equation of state near the phase transition thus focusses the gen- 
eration of anisotropic flow to even earlier times, when the system is still 
entirely partonic and has not even begun to hadronizc. At RHIC energies 
this means that almost all of the finally observed elliptic flow is created dur- 
ing the first 3-4 fm/c of the collision and reflects the hard QGP equation 
of state of an ideal gas of massless particles (c 2 s = | ) 01 Microscopic kinetic 
studies of the evolution of elliptic flow lead to similar estimates for this 

time sca i e raniia 

We close this Section with a beautiful example of elliptic flow from out- 
side the field of heavy-ion physics where the hydrodynamically predicted 
spatial expansion pattern show n in Figure |S] has for the first time been 
directly observed experimentally!^ Figure 1101 shows absorption images of 
an ensemble of about 200,000 6 Li atoms which were captured and cooled 
to ultralow temperatures in a C0 2 laser trap and then suddenly released 
by turning off the laser. The trap is highly anisotropic, creating a pencil- 
like initial spatial distribution with an aspect ratio of about 29 between 
the length and diameter of the pencil. The interaction strength among the 
fcrmionic atoms can be tuned with an external magnetic field by exploiting 
a Feshbach resonance. The pictures shown in Figure correspond to the 
case of very strong interactions. The right panels in Figure ED show that 
the fermion gas expands in the initially short ( "transverse" ) direction much 
more rapidly than along the axis of the pencil. As argued in the paperpSI 
the measured expansion rat es in either direction are consistent with hy- 
drodynamic calculations At late times the gas evolves into a pancake 
oriented perpendicular to the pencil axis. The aspect ratio passes through 
1 (i.e. e x = 0) about 600 /us after release and continues to follow the hydrody- 
namic predictions to about 800 /zs after release. At later times it continues 
to grow, but more slowly than predicted by hydrodynamics, perhaps indi- 
cating a gradual breakdown of local thermal equilibrium due to increasing 
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Fig. 10. Left: False color absorption 
images of a strongly interacting degen- 
erate Fermi gas of ultracold 6 Li atoms 
as a function of time after release from 
a laser trap. Right: Atomic density 
distributions in the initially shorter 
(top) and longer (bottom) directions 
at times 0.4 ms (red, narrowest), 1.0ms 
(blue) and 2.0 ms (green, widest) af- 
ter release from the trap. Reprin ted 
with permission from O'Hara et ail23 
© 2002 AAAS. 



dilution. The authors of the papeJ^I argue that, although the scattering 
among the fermions is very strong by design, it does not seem to be enough 
to ensure rapid local thermalization, and that sufficiently fast healing of 
deviations from local equilibrium caused by the collective expansion might 
require the fermions to be in a superfluid state. 
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4. Experimental observables 

Unfortunately, the small size and short lifetime prohibits a similar direct ob- 
servation of the spatial evolution of the fireball in heavy-ion collisions. Only 
the momenta of the emitted particles are directly experimentally accessible, 
and spatial information must be extracted somewhat indirectly using mo- 
mentum correlations. In the present Section we discuss measurements from 
Au+Au collisions at RHIC and compare them with hydrodynamic calcula- 
tions. Most of the available published data stem from the first RHIC run at 
•y/s = 130 A GeV, but a few selected preliminary data from the 200AGeV 
run in the second year will also be studied. Occasionally comparison will 
be made with SPS data from Pb+Pb collisions at y/s — 17 A GeV. 

This Section is subdivided into two major parts: In Section 14. II we dis- 
cuss single-particle momentum spectra, first averaged over the azimuthal 
emission angle and later analyzed for their anisotropies around the reaction 
plane. These data provide a complete characterization of the momentum- 
space structure of the fireball at freeze-out. In particular the analysis of 
momentum anisotropies yields a strong argument for rapid thermalization 
in heavy-ion collisions and for the creation of a quark-gluon plasma. 

In the second part, Section l4~2l we discuss the extraction of coordinate- 
space information about the fireball at freeze-out from Bose-Einstein in- 
terferometry which exploits quantum statistical two-particle momentum 
correlations between pairs of identical bosons. The general framework of 
this metho d is discussed in the accompanying article by Tomasik and 
WiedemannPSI Here we will discuss specific aspects of Bose-Einstein corre- 
lations from hydrodynamic calculations, in particular their dependence on 
the azimuthal emission angle relative to the reaction plane and its implica- 
tions for the degree of spatial deformation of the fireball at freeze-out. 

4.1. Momentum space observables 

The primary observables in heavy-ion collisions are the triple-differential 
momentum distributions of identified hadrons i as a function of collision 
centrality (impact parameter b): 



We have expanded the dependence on the azimuthal emission angle w„ 
relative to the reaction plane into a Fourier series! 2 ^ Due to reflection 
symmetry with respect to the reaction plane, only cosine terms appear 




(23) 
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in the expansion. As explained before, we restrict our attention to midra- 
pidity, y — hi[(E+p z )/(E— p z )\ = 0, where all odd harmonics (in particular 
the directed flow coefficient v\) vanish. Accordingly, we have also dropped 
the y-dependence of the elliptic flow coefficient v\ but kept its depen- 
dence on t ran sverse momentum and impact parameter. Hydrodynamic 
calculationJSl predict the next higher order coefficient v\ to be very small 
(< 0.1%), and up to now it has not been measured at RHIC. The spectra 
and flow coefficients depend on the hadron species i via the rest mass rrii, 
and this will be seen to play a crucial role. 

As already mentioned, the parameters of the hydrodynamic model 
are fixed by reproducing the measured centrality dependence of the total 
charged multiplicity dN^/dy as well as the shape of the pion and pro- 
ton spectra in central collisions at midrapidity (see below). The shapes of 
other hadron spectra, their centrality dependence and the dependence of 
the elliptic flow coefficient v\ on pj>, centrality and hadron species i are 

|6-| I 

then all parameter free predictions of the mod el f^ 1 The same holds for 
all two-particle momentum correlations! 13 ! 14 ' 74 ! These predictions will be 
compared with experiment and used to test the hydrodynamic approach 
and to extract physical information from its successes and failures. 



4.1.1. Single particle spectra 

The free parameters of the hydrodynamic model are the starting (thermal- 
ization) time r equ , the entropy and net baryon density in the center of the 
reaction zone at this time, and the freeze-out energy density ed oc - The cor- 
responding quantities at other fireball points at r oqu are then determined by 
the Glauber profiles discussed in Sec. 12.41 fsee discussion below Fig. EJ - The 
ratio of net baryon to entropy density is fixed by the measured proton/pion 
ratio. Since the measured chemical composition of the final state at RHIC 
was founG^ to accurately reflect a hadron resonance gas in chemical equi- 
librium at the hadronization phase transition, we require the hydrodynamic 
model to reproduce this p/ir ratio on a hypersurface of temperature T cr it- 
By entropy conservation, the final total charged multiplicity dN^/dy fixes 
the initial product (s • r) e q U P^^^^ The value of r oqu controls how much 
transverse flow can be generated until freeze-out. Since the thermal motion 
and radial flow affect light and heavy particles differently at low p T K>7 l &j l 
a simultaneous fit of the final pion and proton spectra separates the radial 
flow from the thermal component. The final flow strength then fixes r equ 
whereas the freeze-out temperature determines the energy density edec at 
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decoupling. 

The top left panel of Fig.^Jshows the hydrodynamic fiP-^to the trans- 
verse momentum spectra of positive pions and antiprotons, as measured by 
the PHENIX and STAR coll aborations in central (b = 0) Au+Au collisions 
at ^/s = 130 A GeV ! 84 l 85 * 8 -Sl The fit yields an initial central entropy den- 
sity s equ = 95 fm~ 3 at an equilibration time r oqu = 0.6 fm. This corresponds 
to an initial temperature of T equ = 340 MeV and an initial energy density 
e = 25GeV/fm 3 in the fireball center. (Note that these parameters satisfy 
the "uncertainty relation" r cqu • T cqu al.) Freeze-out was implemented on 
a hypersurface of constant energy density with ed cc = 0.075 GeV/fm 3 . 





SPS 


RHIC1 


RHIC2 


y^SNN (GeV) 


17 


130 


200 


Sequ (fm -3 ) 


43 


95 


110 


Tequ (MeV) 


257 


340 


360 


Tequ (fm/c) 


0.8 


0.6 


0.6 



Table 1. Initial conditions for SPS and RHIC energies used to fit the particle spectra 
from central Pb+Pb or Au+Au collisions. s e qu and T equ refer to the maximum values 
at r equ in the fireball center. 

The fit in the top left panel of Fig. 1111 was performed with a chemical 
equilibrium equation of state. Use of such an equation of state implicitly 
assumes that even below the hadronization temperature T cvli chemical equi- 
librium among the different hadron species can be maintained all the way 
down to kinetic freeze-out. With such an equation of state the decoupling 
energy edec = 0.075 GeV/fm 3 translates into a kinetic freeze-out temper- 
ature of Td cc ~130MeV. The data, on the other hand, shcwl^ that the 
hadron abundances freeze out at Tchcm ~ Tcrit, i.e. already when hadrons 
first coalesce from the expanding quark-gluon soup the inelastic processes 
which could transform different hadron species into each other are too slow 
to keep up with the expansion. The measured p/it ratio thus does not agree 
with the one computed from the chemical equilibrium equation of state at 
the kinetic freeze-out temperature Td oc = 130 MeV, and the latter must be 
rescaled by hand if one wants to reproduce not only the shape, but also the 
correct normalization of the measured spectra in Fig. 1111 

A better procedure would be to use a chemical non-equilibrium equation 
of state for the hadronic in which for temperatures T below 

Tchcm ^ ne chemical potentials for each hadronic species are readjusted in 
such a way that their total abundances (after decay of unstable resonances) 
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are kept constant at the observed values. This approach has recently been 
applieo^ to newer RHIC data at */s = 200 A GeV and will be discussed 
below. 
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Fig. 11. Identified pion, antiproton and kaon spectra for 



0.5 



1.5 



p T (GeV) 



PHENIX 



84 91 



and STAR' 



85 86 



7 = 130 GeV from the 
collaborations in comparison with results from a hydro- 
dynamic calculation^! The top left panel shows pion and (anti-)proton spectra from 
central collisions. Shown in the other panels are spectra of five different centralities: 
from most central (top) to the most peripheral (bottom). The spectra are successively 
scaled by a factor 0.1 for clarity. 



For a single hadron species, the shape of the transverse momentum spec- 
trum allows combinations of temperature and radial flow which are strongly 
anticorrelated. 1 ^ By using two hadron species with significantly different 
masses this anticorrelation is strongly reduced albeit not completely elim- 
inated. Consequently, the above procedure still leaves open a small range 
of possible variations for the extracted initial and final parameters. Within 
this range, we selected a value for r oqu which is, if anything, on the large 
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side; some of the hadron spectra would be fit slightly better with even 
smaller values for t qu or a non-zero transverse flow velocity already at T C q U 
(see below). Table 1 su mmariz es the i nitial condit ions applied in our hydro- 
dynamic studies at SPS RfflC l' 13 ' 23 ' 4 ^ and RHIC^ energies. 

Once the parameters have been fixed in central collisions, spectra at 
other centralities and for different hadron species can be predicted without 
introducing additional parameters. The remaining three panels of Fig. 1111 
show the transverse momentum spectra of pions, kaons and antipro tons 
in fiv e diffe rent centrality bins as observed by the PHEN1 j j^4 | 91 | an( j 

STAR MI 

collaborations. For all centrality classes, except the most pe- 
ripheral one, the hydrodynamic predictions (solid lines) agree pretty well 
with the data. The kaon spectra are reproduced almost perfectly, but for pi- 
ons the model consistently underpredicts the data at low p^. This has now 
been understood to be largely an artifact of having employed in these calcu- 
lations a chemical equilibrium equation of state all the way down to kinetic 
freeze-out. More recent calculations^ with a chemical non-equilibrium 
equation of state, to be compared to s/s — 200 A GeV data below, show 
that, as the system cools below the chemical freeze-out point T c hem ~ T CI a , 
a significant positive pion chemical potential builds up, emphasizing the 
concave curvature of the spectrum from Bose effects and increasing the 
feeddown corrections from heavier resonances at low pt- The inclusion of 
non-equilibrium baryon chemical potentials to avoid baryon-antibaryon an- 
nihilation further amplifies the resonance feeddown for pions. 

Significant discrepancies are also seen at large impact parameters and 
large transverse momenta pt > 2.5 GeV/c. This is not surprising since high- 
Pt particles require more rescatterings to thermalize and escape from the 
fireball before doing so. This is in particular true in more peripheral colli- 
sions where the reaction zone is smaller. 

For the calculations shown in Fig.^Jthe same value edec was used for all 
impact parameters. Phenomenological studies^ using a hydrodynamically 
motivated parametrizatiorP^ to describe pion and antiproton spectra from 
200 A GeV Au+ Au collisions in a large number of centrality bins indicate 
somewhat earlier freeze-out, at higher temperature and with less transverse 
flow, in the most peripheral collisions (see Fig. I12|) . In the hydrodynamic 
model this can be accommodated by allowing the freeze-out energy density 
edcc to increase with impact parameter. A consistent determination of edec 
from the kinetic decoupling criterion is expected to automatically yield such 
a behavior. Such a calculation would use the fitted value for edcc extracted 
from central collision data to determine the unknown proportionality con- 
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stant between the local expansion and scattering time scales at decoupling 
(see discussion in Sees. 12.51 and lH.l} . and then calculate edec for other impact 
parameters by using the kinetic freeze-out criterion T exp = n r scat t with the 
same constant n extracted from central collisions. So far this has not been 
done, though. 



PHENIX Preliminary 




Fig. 12. Kinetic freeze- 
out temperature Tf G = Tdec 
and transverse flow veloc- 
ity at the fireball edge, 
extracted from a simulta- 
neous fit of a flo w sp ectrum 
parametrizatiorlSJJ to 7T , 
if , p and p spectra from 
200AGeV Au+Au colli- 
sions over the entire ra nge 
of collision centralitiesl22l 
More peripheral collisions 
(small numbers -/V par t of 
participating nucleons) are 
seen to decouple earlier, at 
higher freeze-out tempera- 
ture and with less trans- 
verse flow. 
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Without transverse flow, thermal spectra exhibit rriT-scalin^^, i.e. af- 
ter appropriate rescaling of the yields all spectra collapse onto a single 
curve. Transverse collective flow breaks this scaling at low px < m o (i- e - f° r 
non-relativistic transverse particle velocities) by an amount which increases 
with the particle rest mass mo When plotting the spectra against 

Pt instead of mx, any breaking of mx-scaling is at least partially masked 
by a kinematic effect at low px which unfortunately again increases with 
the rest mass mo- To visualize the effects of transverse flow on the spec- 
tral shape thus requires plotting the spectrum logarithmically as a func- 
tion of mx or mx— m n- Such plots can be found in recent experimental 
, 94 95 96 , 



publications, 



and although the viewer's eye is often misled by su- 



perimposed straight exponential lines 



-m T /T B{f 



, a second glance shows 



a clear tendency of the heavier hadron spectra to curve and to begin to 
develop a shoulder at low transverse kinetic energy mx— mp, as expected 
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from transverse flow. 

One such example is shown in Fig. 1131 where preliminary spectra of il 
hyperonJ^ are compared with hydrodynamic predictions. For this com- 
parison the original calculations for 130 A GeV Au+Au collisions^ were 
repeated with RHIC2 initial conditions and a chemical non-equilibrium 
equation of state in the hadronic phaseP^I The solid lines are based on 
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Fig. 13. Transve rse mass spectrum of Q hyperons from central 200 A GeV Au+Au 
collisions at RHICl^l The curves are hydrodynamic calculations with different initial 
and freeze-out conditions: Solid lines correspond to the default of no initial transverse 
flow at T cqu , dashed lines assume a small but non-zero radial flow, v r = tanh(ar) with 
a = 0.02 fm -1 , already at T cqu . The lower (thin) set of curves assumes fi-decoupling at 
T" cr it = 164 MeV, the upper (thi ck) set of curves decouples the Q together with the pions 
and protons at T dcc = 100 MevPSl 

default parameters (see Table 1) without any initial transverse flow at r oqu . 
(The dashed lines will be discussed further below.) Following a suggestion 
that Q hyperons, being heavy and not having any known strong coupling 
resonances with pions, should not be able to participate in any incre ase 
of the radial flow during the hadronic phase and thus decouple earlyj*^ 
we show two solid lines, the steeper one corresponding to decoupling at 
edcc = 0.45 GeV/fm 3 , i.e. directly after hadronization at T cr j t , whereas the 
flatter one assumes decoupling together with pions and other hadrons at 
edcc = 0.075 GeV/fm 3 . The data clearly favor the flatter curve, suggesting 
intense rescattering of the f2's in the hadronic phase. The microscopic mech- 
anism for this rescattering is still unclear. However, without hadronic rescat- 
tering the hydrodynamic model, in spite of its perfect local thermalization 
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during the early expansion stages, is unable to generate enough transverse 
flow to flatten the f2 spectra as much as required by the data. Partonic 
hydrodynamic flow alone can not explain the spectrum. 

We close this subsection with a comparison of pion, kaon and hadron 
spectra from 200 A GeV Au+Au collisions (RHIC2) with recent hydrody- 
namic calculations which correctly implement chemical decoupling at T cr ; t . 
Figure sho ws a compilatio n of preliminary spectra from the four RHIC 
collaborations ! I I I ^For better visibility, the 7r _ , K~ and p spectra 
are separated artificially by scaling factors of 10 and 100, respectively. The 
lines reflect hydrodynamic results. In these calculations the particle num- 



PHENIXi 
STAR [ 
PHOBOS | 
BRAHMS I 




P T (GeV) 

Fig. 14. Particle spectra of tt~ , K~ an d antiprotons at -^/snn = 200 GeV as measured 
by the four large experiments at RHIcMIMEOl The lines show hydrodynamic 
results under various considerations (see text) ESI 



bers of all stable hadron species were conserved throughout the hadronic 
resonance gas phase of the evolution, by introducing appropriate chemi- 
cal potentials! 9 ' I w l w l It turns out that such a chemical non-equilibrium 
equation of state has almost the same relation p(e) between the pressure 

and energy density as the equilibriu m o ne, and that the hydrodynamical 

IAVIAqI 

evolution remains almost unaltered! 9 ' l qy l However, the relation between the 
decoupling energy density edec and the freeze-out temperature Td ec changes 
significantly, since the non-equilibrium equation of state prohibits the anni- 
hilation of baryon-antibaryon pairs as the temperature drops. Consequently, 
at any given temperature below T c h e m the non-equilibrium equation of state 
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contains more heavy baryons and antibaryons than the equilibrium one and 
thus has a higher energy density. 

The same energy density edcc = 0.075 GeV/fm 3 then translate s into a 
significantly lower freeze-out temperature T& ec w 100 MeV! l °^ l Da l The cor- 
responding results are given as thick solid (red) lines in Fig. El The thin 
solid (blue) lines in the Figure, shown for comparison, were calculated 
by assuming kinetic freeze-out already at hadronization, T cr j t = 165 MeV. 
The data clearly favor the additional radial boost resulting from the 
continued buildup of radial flow in the hadronic phase. Still, even at 
edcc = 0.075 GeV/fm 3 , the spectra are still steeper than the data and the 
previous calculations with a chemical equilibrium equation of state shown 
in Fig. 1111 reflecting the combination of the same flow pattern with a lo wer 
freeze-out temperature. Somewhat unexpectedly, the authors of the study^l 
were unable to significantly improve the situation by reducing edcc even fur- 
ther: The effects of a larger radial flow at lower edcc were almost completely 
compensated by the accompanying lower freeze-out temperature, leading to 
only modest improvements for kaons and protons and almost none for the 
pions. 

This motivated the author to introduce a small but non-vanishing 
transverse "seed" velocity already at the beginning of the hydrodynamic 
stage. The dashed lines in Fig. El (and also earlier in Fig. 113(1 show hy- 
droynamic calculations with an initial transverse flow velocity profile given 
by v r (r, r e q U ) = tanh(ctr) with a = 0.02 fm _1 . This initial transverse kick 
is seen to significantly improve the agreement with the pion, kaon and 
antiproton data up to px > 1-5 — 2GeV/c for pions and kaons and up to 
Pt si 3.5GeV/c for (anti)protons~3 It can be motivated by invoking some 
collective (although not ideal hydrodynamic) transverse motion of the fire- 
ball already during the initial thermalization stage, although the magnitude 
of the parameter a requires further study. 

In Figure El the kaon and antiproton spectra were divided by factors 
of 10 and 100, respectively, for better visibility. If this is not done one no- 
tices that the antipr oton s pectrum crosses the pion spectrum at around 
Pt ~ 2 — 2.5 GeV/cj ^^ l i.e. at larger pt antiprotons are more abun- 
dant than pions! This became known as the u p/it~ > 1 anomaly" and has 
attracted significant attention! -*-^ ! Here the word "anomaly" arises from a 
comparison of this ratio in central Au+ Au with pp collisions and with string 
fragmentation models which both givep/7r~ ratios much below 1. However, 
in Au+Au collisions string fragmentation is expected to explain hadron pro- 
duction only at rather large pt! ^^ and in the hydrodynamic picture which 
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is successful at px < 2 GeV/c there is actually nothing anomalous about a 
p/ir ratio that exceeds 1. 

To see this let us first look at a thermal system without flow. The 
corresponding transverse mass spectra are in good approximation simple 
exponentials in %jf- whose ratios at fixed mx are simply given by the ratios 
{§pl 9-K) e ^ >lp ~^ 7 '^ T °f their spin-isospin degeneracies and fugacities. For suf- 
ficiently large px 2> fn p , mx ~Pt and the same holds true for the ratio of 
the px-spectra at fixed px- It still holds true in the presence of transverse 
flow which, at sufficiently large px, simply flattens all mx-slopes by a com- 
mon blucshift factoJ^ESESl VlSit, Since antiprotons have a 2-fold spin 
degeneracy and pions have none, we see that the asymptotic hydrodynamic 
p/n ratio is above unity if the chemical potentials are sufficiently small. 
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100 MeV 


0.7 
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Table 2. Upper part: Chemical potentials of protons, antiprotons, pions and kaons at 
the chemical (164 MeV) and kinetic (100 MeV) freeze-out temperatures for 200AGCV 
Au+Au collisions EH Lower part: Asymptotic particle ratios for these hadrons at fixed 
large pxi for two assumed hydrodynamic freeze-out temperatures of 164 and 100 MeV, 
respectively (see text). 

In Au+ Au collisions at v / snn — 200 GeV the baryon chemical potential 
at chemical freeze-out is smalP^ (29 MeV) and the pion chemical potential 
vanishes. Correspondingly, the asymptotic p/ir + and p/tt~ ratios are both 
close to 2 (see Table 2). As the system cools below the chemical freeze- 
out temperature, however, pions, kaons and both protons and antiprotons 
all develop significant positive chemical potentials which are necessary to 
keep their total abundances (after unstable resonance decays) fixed at their 
chemical freeze-out values^ (second row in Table 2). As a consequence, 
the asymptotic p/ir and p/n ratios increase dramatically, to 40 and 28, 
respectively, and even the asymptotic K + /ir + ratio increases from 1 to 2.4. 

We see that cooling at constant particle numbers strongly depletes the 
pions at high px in favor of high-px baryons and kaons. Correspondingly, 
the ratios of the hydrodynamic spectra, shown in Fig. rise far above 
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Fig. 15. Hydrodynamic predictions for the midrapidity ratios of p/ir + , p/ir~, and 
J\+/-7r+ as functions of px f° r 200 A GeV Au+Au collisions, extracted from the top 
pairs of curves shown in Fig. 1141 (see description in text). 

unity at large pr- We should stress that the increase with px of these ratios 
at small px is generic for thermalized spectra and independent of whether 
or not there is radial flow. It is a simple kinematic consequence of replotting 
two approximately parallel exponentials (more exactly: K2-functions) in mx 
as functions of pt and taking the ratio. Due to the larger rest mass the px- 
spectrum of the heavier particle is flattened more strongly at low px than 
that of the lighter particle, yielding for their ratio a rising function of px- 
The additional flattening from radial flow, which again affects the heavier 
particles more strongly than the light ones, further increases this tendency. 

It is worth pointing out that a p/ir ratio well above 2 and and a K/tt 
ratio well above 1, as hydrodynamically predicted for px >2.5GeV/c (see 
Fig. I15fl would, when taken together with the measured global thermal 
yields, provide a unique proof for chemical and kinetic decoupling happen- 
ing at different temperatures. Unfortunately, as we will see in more detail 
later, hydrodynamics begins to seriously break down exa ctly in this inter- 
esting pt domain, and the experimentally observed ratioJ^ never appear 
to grow much beyond unity before decreasing again at even higher px , even- 
tually perhaps app roaching the small asymptotic values expected from jet 
fragmentation^^ 

4.1.2. Mean transverse momentum and transverse energy 

The good agreement of the hydrodynamic calculations with the experimen- 
tal transverse momentum spectra is reflected in a similarly good description 
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of the measured average transverse momenta!^ Figure El snows a compar- 
ison of (pt) for identified pions, kaons, pro tons an d antiprotons measured 
by PHENIX in 200 A GeV Au+Au collisionPdMil with the hydrodynamic 
results.ESl The bands reflect the theoretical variation resulting from pos- 
sible initial transverse flow already at the beginning of the hydrodynamic 
expansion stage, as discussed at the end of the previous subsection. The 
figure shows some discrepancies between hydrodynamics and the data for 
peripheral collisions (small -/Vp art ) which are strongest for the kaons whose 
spectra are flatter at large impact parameters than predicted by the model. 
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Fig. 16. Mean tr ansvers e momenta of pions, kaons and (anti)protons from 200 A GeV 
Au+Au collisions | 99 | 1Q3 Hydrodynamic results are included as bands whose lower ends 
reflect an initialization without initial transverse flow while the up per ends correspond 
to an initial transverse flow field v r = tanh(ar) with a = 0.02 fm — 

FigureElshows the total transverse energy per emitted charged hadron 
as a functi on o f collision centrality. The data are from Pb+ Pb collis ions 
at the SP922D an d Au+Au collisions at two RH1C energiesPSIIlilSI Al- 
though both the charged particle multiplicity and total transverse energy 
vary strongly with the number of participating nucleons (and one or the 
other are therefore often used to determine the collision centrality), the 
transverse energy per particle is essentially independent of the centrality. It 
also depends only weakly on the collision energy. 

The superimposed band in Figure lTTI reflects hydrodynamic calculations 
for Au+Au collisions at ^fs = 200 A GeV with and without initial transverse 
flow, as before. The slight rise of the theoretical curves with increasing Ap art 
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Fig. 17. Transverse energy per charged hadron as a functio n of collision centrality, for 
Au+Au and Pb+Pb collisions at th ree di fferent beam energies l^Wl+W l lvv lSupcrimposcd 
on the original expe rim ental FigurJi^SJ are hydrodynamic results for Au+Au collisions 



200 ( icY."-' The lower end of the band results from an initialization without 



initial transverse flow, the upper end reflects an initial transverse flow field v r = tanh(or) 
with a = 0.02 fin -1 . 

can be attributed to the larger average transverse flow developing in more 
central collisions, resulting from the higher initial energy d ensi ty and the 
somewhat longer duration of the expansion until freeze-out!^ Successful 
reproduction of the data requires a correct treatment of the chemical com- 
position at freeze-out (by using a chemical non-equilibrium hadron equa- 
tion of state below T cr ; t ). If one instead assumes chemical equilibrium of 
the hadron resonance gas down to kinetic freeze-out, hydro dyna mics over- 
predicts the transverse energy per particle by about 15-20%!^ 



4.1.3. Momentum anisotropics as early fireball signatures 

In non-central nuclear collisions, or if the colliding nuclei are deformed, the 
nuclear overlap region is initially spatially deformed (see Fig. |3J). Interac- 
tions among the constituents of the matter formed in that zone transfer this 
spatial deformation onto momentum space. Even if the fireball matter does 
not interact strongly enough to reach and maintain almost instantaneous 
local equilibrium, and a hydrodynamic description therefore fails, any kind 
of re-interaction among the fireball constituents will still be sensitive to the 
anisotropic density gradients in the reaction zone and thus redirect the mo- 
mentum flow preferably into the direction of the strongest density gradients 
(i.e. in the "short" direction) | 75 | 76 | 78 | 1(J9| resu i-t; j s a momentum-space 
anisotropy, with more momentum flowing into the reaction plane than out 



February 5, 2008 7:22 



WSPC/Trim Size: 9in x 6in for Review Volume 



42 Hydrodynamic description of ultrarelativistic heavy-ion collisions 

Of it. 

Such a "momentum-space reflection" of the initial spatial deformation 
is a unique signature for re-interactions in the fireball and, when observed, 
proves that the fireball matter has undergone significant nontrivial dynam- 
ics between creation and freeze-out. Without rescattering, the only other 
mechanism with the ability to map a spatial deformation onto momentum 
space is the quantum mechanical uncertainty relation. For matter confined 
to smaller spatial dimensions in x than in y direction it predicts Ap x > Ap y 
for the corresponding widths of the momentum distribution. However, any 
momentum anisotropy resulting from this mechanism is restricted to mo- 
menta p ~ l/(size of the overlap zone) which for a typical fireball radius of a 
few fm translates into a fraction of 200 MeV/ c. This is the likely mechanism 
for the momentum anisotropy observed I in calculations of the classical 
dynamical evolution of a postulated deformed "color glass condensate" cre- 
ated initially in the collision. Unlike the experiment al dat a, this momentum 
anisotropy is concentrated around relatively low pt .^23 

Whatever the detailed mechanism responsible for the observed momen- 
tum anisotropy, the induced faster motion into the reaction plane than 
perpendicular to it ( "elliptic flow" ) rapidly degrades the initial spatial de- 
formation of the matter distribution and thus eliminates the driving force 
for any further increas e of the anisotropic flow. Elliptic flow is therefore 
"self-quenching" 1 75 1 76 1 and any flow anisotropy measured in the final state 
must have been generated early when the collision fireball was still spatially 
deformed (see Fig. EJ) . If elliptic flow does not develop early, it never devel- 
ops at all (see also Sec. 14.1 3} . It thus reflects the pre ssure and stiff ness of 
the equation of state during the earliest collision stages E 5 l 7t? l 77l4l Z51but ( m 
contrast to many other early fireball signatures) it can be easily measured 
with high statistical accuracy since it affects all final state particles. 

Microscopic kinetic models show (see Fig. 1191 below ) that, for a given 
initial spatial deformation, the induced momentum space anisotropy is a 
monotonically rising function of the strength of the interaction among the 
matter constituents P^^^ The maximum effect should thus be exp ected 
if their mean free path approaches zero, i.e. in the hydrodynamic limit l^iiSl 
Within this limit, we will see that the magnitude of the effect shows some 
sensitivity to the nuclear equation of state in the early collision stage, but 
the variation is not very large. This implies that, since the initial spatial 
deformation can be computed from the collision geometry (the average 
impact parameter can be determined, say, from the ratio of the observed 
multiplicity in the event to the maximum multiplicity from all events), the 
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observed magnitudes of the mom entum an isotropics, and in particular their 
dependence on collision centrahtyP^SIIIIl provide valuable measures for the 
degree of thermalization reached early in the collision. 

Experimentally this program was first pursued at the SPS in 158 A GeV 
Pb+Pb collisions! ^ ^ These data still showed significant sensitivity 
to details of the analysis procedure^EI anc j thus remained somewhat 
inconclusive!^! Qualitatively, the SPS data (where the directed and elliptic 
flow coefficients, V\ and 1)2, can both be measured) confirmed Ollitrault's 
1992 prediction'" that near midrapidity the preferred flow direction is into 
the reaction plane, supporting the conclusions from earlier measurements 
in Au+Au collisions at the AGS ^^ I where a transition from out-of-plane 
to in-plane elliptic flow had been found between 4 and 6 A GeV beam en- 
ergy. A comprehensive quantitative discussion of elliptic flow became first 
possible with RHIC data, because of their better statistics and improved 
event plane resoluti on (d ue to the larger event multiplicities) and also as 
a result of improve d I analysis techniques. a In the meantime the latter 
have also been re-applied to SPS data and prod uced very detailed results 
from Pb+Pb collisions at this lower beam energy! ' ™ 

4.1 .4. Elliptic flow at RHIC 

The second published and still among the most important results from 
Au+Au collisions at RHIC was the centrality and pt dependence of the 
elliptic flow coefficient at midrapidity! - ^ For central to midperipheral col- 
lisions and for transverse momenta pt^1-5 GeV/c the dat a, wer e found 
to be in stunning agreement with hydrodynamic predictions j ^ l ^ l as seen 
in Fig. ^1 In the left panel, the ratio n c h/n max of the charged particle 
multiplicity to the maximum observed value is used to characterize the 

collision centrality, with the most central collisions towards the right near 

[Toil 

1. n c h/n max = 0.45 corresponds to an impact parameter bfa 7fm! *Up to 
this value the observed elliptic flow V2 is found to track very well the in- 
creasing initial spatial deformation e x of the nuclear overlap zone|l^ as 



a Onc of the important experimental issues is the separatiorl ^ ^ 1 1 S3 of collective "flow" 
contributions to the observed momentum anisotropics from "non-flow" angular cor- 
relations, su ch as Bo se-Einstein correlations )-^ ^ 1 correlations arising from momentum 
conservation 

in Kin »i 

and two-particle momentum correlations from resonance decays 
and jet productionl ^^^ l Many of the data shown in this r eview have not yet been cor- 
rected for non-flow contributions, but subsequent analysiJ ^^^ I has shown that for the 
RHIC data and in the p-r-range of interest for this review these corrections are small 
(<15%). 
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predicted by hydrodynamics.^ 




Fig. 18. Elliptic flow of unidentified charged particles in 130 AGeV Au+Au collisions, 
integrated over px as function of centrality (left) and for minimum bias collisions as a 
function of pt (right). Both data sets (symbols with e rror bars) are from the original 
STAR publicatiomESl The vertical bars in the left pane(i^5] indicate the range of earlier 
hydrodynamic predictions for a v arie ty of equations of state and initial conditions0The 
top three curves in the right panel^l represent hydrodyna mic p redictions for semiperiph- 
eral collisions with initial conditions tuned to the observed^! total charged multiplicity 
in central collisions where V2 vanishe s. D ifferent curves correspond to different equations 
of state and freeze-out temperatures!^! 



Following this discovery it soon became clear that the agreement of the 
data with the hydrodynamic calculations could not be accidental and in fact 
allows to draw a number of very strong and important conclusions. These 
conclusions refer to soft particle production, that is to mesons with trans- 
verse momenta up to about 1.5GeV/c and baryons with <2.5GeV/c. 
This momentum range covers well over 99% of the produced particles. This 
means that we are talking about the global dynamical features of the bulk of 
the fireball matter. Of course, the small fraction of particles emitted with 
larger transverse momenta carry very important information themselves, 
but their behavior is not expected to be controlled by hydrodynamics in 
the first place, and they are not the subject of our discussion here. It should 
be noted, however, that interpreting the behavior of high-pT particles does 
require a prior understanding of the global fireball dynamics which is the 
subject of this review. In the following we discuss the aforementioned con- 
clusions, as well as a number of additional theoretical and experimental 
aspects. 
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Strong rescattering: 

It was quickly realized^ that the measure almost linear rise of 
the charged particle (i.e. predominantly pionic) elliptic flow with px re- 
quires strong rescattering among the fireball constituents. Figure ITP1 shows 
the results from microscopic simulations which describe the dynamics of 
the early expan sion stage by solving a Boltzmann equation for colliding 
on-shell partonsP^ The different curves are parametrized by the transport 
opacity £ = uodNg/dn involving the product of the parton rapidity density 
and cross section in the early collision stage. As the opacity is increased, 
the elliptic flow is seen to approach the data (and the hydrodynamic limit) 
monotonically from below. Whereas the hydrodynamic limit predicts a con- 
tinuous rise of V2(pt), the elliptic flow from the parton cascade saturates 
at high pt, as also seen in the datcP25l ( see Fig. below). This is due to 
incomplete equilibration at high p^: the critical pt at which the cascade 
results cease to follow the hydrodynamic rise shifts to higher (lower) values 
as the transport opacity is increased (decreased). 
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Fig. 19. Impact parameter averaged elliptic flo w as a function of transverse momen- 
tum. The exper ime ntal data points from STAr4i^l are compared with parton cas- 
cade calculationJZSl with varying transport opacities £. Figure taken from Molnar and 
GyulassyEHl 



It is interesting to observe that stronger rescattering manifests itself 



February 5, 2008 7:22 



WSPC/Trim Size: 9in x 6in for Review Volume 



46 Hydrodynamic description of ultrarelativistic heavy-ion collisions 

in this way, i.e. by following the hydrodynamic curve with the full slope to 
higher pt and not by producing a hydro-like quasi-linear rise with a reduced 
slope. In view of this the elliptic flow data at large im pact p arameters (see 
Fig. E3 below) and at lower collision energieJ^I 11 ^ 122 ! 1 ^ , which show a 
linear rise of V2(pt) with a smaller slope than hydrodynamically predicted, 
pose an unresolved puzzle which is not simply explained by incomplete local 
thermalization. 

Figure IT§1 also shows that very high transport opacities are necessary if 
the parton cascade is required to follow the data to pt ~ 1.5—2 GeV/c. The 
necessary values exceed perturbative expectations by at least a factor 30,E3 
raising the question which microscopic interactio n mechanism is responsi- 
ble for th e lar ge observed elliptic flow!E22II2EI However, it was recently 
discovered 'that the partonic elliptic flow may not necessarily have to fol- 
low the hydrodynamic prediction all the way out to pt — 1-5 GeV/c: If the 
elliptic flow of the partons gets transferred to the hadrons by a momentum- 
space coalescence mechanism PI it is sufficient if it behaves hydrodynami- 
cally up to pt ~ 0.7—0.8 GeV/c for the pion and proton v 2 to "look hydro- 
dynamic" up to pr~1.5GeV/c and 2.2-2.4 GeV/c, respectively!^ This 
takes aw ay so me of the pressure for anomalously large partonic transport 
opacities II2xl 

Centrality dependence of elliptic How: 

Figure IT51 showed that in peripheral collisions the pr-urtegrated elliptic 
flow lags behind the hydrodynamic predictions. This may reflect incomplete 
thermalization in the smaller fireballs created in these cases. To study this 
in more detail, Fig. 1201 shows th e PT -differential elliptic flow for pions and 
protons for three centrality bins! Due to limitations for particle identi- 
fication, the data cover only the low-px region up to about 800MeV/c. In 
this pt region, the left panel shows good agreement of V2{pt) for pions with 
the hydrodynamic predictionJ^ for central and midcentral collisions, but 
smaller elliptic flow than predicted for the most peripheral bin (45-85 % cen - 
trality, corresponding to an average impact parameter of about 11 fm^l). 
The graph clarifies that for peripheral collisions the smaller-than-predicted 
PT-integrated elliptic flow seen in Fig. 1181 arises mostly from a smaller-than- 
predicted slope of the pr- differential elliptic flow for pions. For the most 
peripheral collisions this slope is about 20% less than expected if also there 
the reaction zone were able to fully thermalize. In view of the large aver- 
age impact parameter in this centrality bin, it is rather surprising that the 
discrepancy to the hydrodynamic predictions is not larger. 
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Fig. 20. Elliptic flow of pions (left) and (anti-)protons (right) as measured by the STAR 
CollaboratioiffiS D in three different centrality bins, fncluded are results from a hydrody- 
namic predictionl^U 



The right panel in Fig. 1201 shows a similar comparison for protons and 
antiprotons. Due to the limited statistics of the data, which also were not 
fully corrected for feeddown from weak A decays, no strong conclusions 
can be drawn from the plot, but the data seem to be generally on the low 
side compared to the hydrodynamic curves. However, this can have other 
reasons than a breakdown of hydrodynamics, due to a specific sensitivity 
of the elliptic flow of heavy hadrons to the nuclear equation of state (see 
below) . 
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Fig. 21. Left: The px _m tegrated elliptic flow V2 scaled by the initial spatial eccentric- 
ity e x as a function of impact parameter. Shown are hydrodynamic results with initial 
conditions corresponding to AGS, SPS and RHIC energies. Quoted are the respective 
values eo for the central energy density in 6 = Pb +P b collisions: eo = 4. 5 Ge V/fm 3 
(AGS), 9.0 GeV/fm 3 (SPS) and 25 GeV/fm 3 (RHIC) Right: PHENIX datJSSIfor the 
ratio A2 = V2/e x measured in 130 A GeV Au+Au collisions at RHIC in two different pt 
regions. Central collisions are to wards the left, peripheral collisions towards the right of 
the diagram. (Figure taken from ESI). 
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In hydrodynamic calculations the finally observed elliptic flow is essen- 
tially proportional to the initial spatial eccentricity e x of the reaction zone 
f Section |2. 4(1 . This is displayed in the left panel of Fig. 1211 where the ellip- 
tic flow scaled by the initial eccentricity is plotted as a function of impact 
paramete3 (note the suppressed zero on the vertical axis). The slight vari- 
ation of the ratio A2 —vijt x with impact parameter reflects changes in the 
stiffness of the equation of state, resulting from the quark-hadron phase 
transition, which are probed as the impact parameter (and thus the initial 
energy density in the center of the reaction zone) is varied^ This will be 
discussed in more detail when we describe the beam energy dependence of 
elliptic flow. 

The righ t pa nel of Fig. |2T] shows RHIC data from the PHENIX 
CollaboratiorPSni for the same ratio, at low and high transverse momenta. 
While the low-pT data agree with the hydrodynamic prediction of an ap- 
proximately constant ratio V2 /e x , at high pt the scaled elliptic flow is seen 
to decrease for more peripheral collisions. This is consistent with the earlier 
discussion of a gradual breakdown of hydrodynamics for increasing pt and 
impact parameter. 

Elliptic How of different hadron species: 

Hydrodynamics predicts a clear mass-ordering of elliptic flowPl As the 
collective radial motion boosts particles to higher average velocities, heavier 
particles gain more momentum than lighter one s, le ading to a flattening of 
their spectra at low transverse kinetic energies!""' When plotted against 
Pt this effect is further enhanced by a kinematic factor arising from the 
transformation from tot to pt (see Fig. El an d discussion below Fig. | 12H . 
This flattening reduces the momentum anisotropy coefficient V2 at low j?t]^ 
and the heavier the particle the more the rise of U2(pt) is shifted towards 
larger pt (see top left panel in Fig. ■ 

This effect, which is a consequence of both the thermal shape of the 
single-particle spectra at low pt and the superimposed collective radial flow, 
has been nicely confirmed b y the experim ents: Figure|221 an d the right panel 
of Fig.EHshow that the follow the predicted mass ordering 

out to transverse momenta of about 1.5GeV/c. For K® and A + A much 
more accurate data than those shown in Fig. 1221 ft op rig ht) have recently be- 
come available from 200 A GeV Au+Au collisions) ™ again in quantitative 
agreement with hydrodynamic calculations up to pt — 1.5 GeV/c for kaons 
and up to pt — 2.5 GeV/c for A + A. The inversion of the mass-ordering in 
the data at large pt is caused by the mesons whose «2(pt) breaks away 
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Fig. 22. The pT-diffcrcntial elliptic flow V2(pt) from minimum bias Au+Au collisions 
at RHIC, for differ ent identified hadron species. Top left: hydrodynamic predictions 
for v / i=130AGevlHII Top right: V2(pt) for charged hadrons (mostly pions), K® and 
A + A at y/s = 130 A GeV measu red by the STAR Collaboration and shown together with 
hydrodynamic calculations^^] Bottom row: preliminary results for i>2(pt) of identified 
pions, kaons and protons with n egative (left) and positive (right) charge, measured by 
the PHENIXCpllaboratiorlSll at y/s = 200 A GeV and compared with hydrodynamic 
calculations 



from the hydrodynamic rise and begins to saturate at pt> 1.5GeV/c In 
contrast, baryons appear to behave hydrodynamically to pt — 2.5 GeV /c, 
breaking away from the flow prediction and saturating at significantly larger 
Pt than the mesons. This is consistent with the idea that the partonic el- 
liptic flow established before hadronization exhibits a hydrodynamic rise at 
low pt followed by saturation above pt ~ 750 — 800 MeV/c, and that these 
features are transferred to the observed hadrons by quark coalescenc e, ma n- 
ifesting themselves there at twice resp. three times larger px-values! ' 
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Sensitivity to the equation of state: 

The experimental determination of the nuclear equation of state at high 
densities relies on detailed studies of collective flow patterns generated in 
relativistic heavy- ion collisions^] Since elliptic flow builds up and saturates 
early in the collision, it is more sensitive to the high density equation of state 
than the azimuthally averaged radial flow. Hydrodynamic calculations 
allow to study in the most direct way the influence of the phase transition 
and its strength (the latent heat Aei a t) on the generated flow patte rns. This 
was investigated systematically and in great detail by Teaney et al. j l^ l l^ EOI 
using hydrodynamics to describe the early quark-gluon plasma expansion 
stage (including the quark-hadron phase transition), followed by a kinetic 
afterburner which simulates the subsequent hadronic evolution and freeze- 
out with the relativistic hadron cascade RQMDpSI 

One of their important 
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Fig. 23. Left: px _m t e g ra t e d elliptic flow for 130 A GeV Au+Au as a function o f col lision 
centrality from a hydrodynamic calculation with a hadron cascade afterburnerBSl LH4, 
LH8 and LH16 label three different phase transitions with latent heats Aei pj- = .4. 0.8 
and 1.6 GeV/fm 3 , respectively. Data are from the STAR Collaboration!^ Right: 
Px-differential elliptic flow of identified pions and (anti)protons from minimum bias 
130 A GeV Au+Au collisio ns! -*-^ ! compared with hydrodynamic calculations using Coo- 
per- Frye freeze-out See text for details. 

results is shown in the left panel of Fig. 1231 which (similar to Fig. Ii8jl gives 
the pT-integrated elliptic flow as a function of the normalized particle yield 
as a centrality measure. The three curves correspond to equations of state 
with a first order quark-hadron transitio n, w ith latent heat values Aei a t 
of 0.4, 0.8 and 1.6GeV/fm 3 , respectively!^ Comparison with the results 
from the STAR Collaboration shows that a phase transition of significant 
strength (Aei at > 1 GeV/fm 3 ) is necessary to reproduce the data. Without 
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the softening of the equation of state induced by the phase transition, the 
single-particle spectra are too fiat and the pT-mtegrated elliptic flow comes 
out too large, even though V2(pt) has roughly the correct slope for pions. 

On the other hand, if one eliminates the "QGP push" entirely and 
replaces the hard quark-gluon plasma equation of state above the phase 
transition by a softer hadron resonance gas without phase transition 
(EOS H) , one un derpredicts the hydrodynamic mass-splitting of the el- 
liptic flow! 14 ! 135 ' This is seen in the right panel of the figure which shows 
the pT-diffcrential elliptic flow for pions and protons both with a realistic 
equation of stat e (E OS Q) and for a pure hadron resonance gas (EOS H). 
The proton datcOH 

shown in this plot obviously favor EOS Q, irrespective 
of moderate variations of the freeze-out temperature, indicated by the three 
lines labelled by EOS Q™H T eaney et alfelcame to similar conclusions; 
this implies that the details of how the hadronic resc attering stage is de- 
scribed (hydrodyna micaUy wit h Cooper- Frye freeze- kinetically 

11 AI1 mI^AI 

via a hadron cascad e ! 1 1 D do not matter. 

Rapid thcrmalization: 

We can summarize the comparison between RHIC data and the hydro- 
dynamic model up to this point by stating that hydrodynamics provides 
a good description of all aspects of the single particle momentum spectra, 
from central and semicentral Au+Au collisions up to impact parameters 
&~ 10 fm, and for transverse momenta up to 1.5GeV/c for mesons and up 
to 2.5GeV/c for baryons. Since this px-range covers well over 95% of the 
emitted particles, it is fair to say that the bulk of the fireball matter cre- 
ated at RHIC behaves hydrodynamicaUy, with little indication for non-ideal 
(viscous) effects. As explained, the successful description of the data by the 
hydrodynamic model requires starting the hydrodynamic evolution no later 
than about lfm/c after nuclear impact. This estimate is even conservative 
since it does not take into account any transverse motion of the created 
fireball matter between the time when the nuclei first collide and when the 
fireball has thermalized and the h ydro dy namic expansion begins. We will 
now give an independent argumenPD^2 why thermalization must happen 
very rapidly in order for the elliptic flow signal to be as strong as observed 
in the experiments. 

As shown earlier in this Section, the hydrodynamicaUy predicted elliptic 
flow is proportional to the initial spatial eccentricity e x (Te qu ) at the begin- 
ning of the hydrodynamic evolution. If thermalization is slow, the matter 
will start to evolve in the transverse directions even before r cqu , following its 
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initial locally isotropic transverse momentum distribution. Even if no rein- 
teractions among the produced particles occur, this radial free-streaming 
motion dilutes the spatial deformation, although not quite as quickly as in 
the opposite limit of complete thermalization where it decreases faster due 
to anisotropic hydrodynamic flow (see Fig.^J). Thus, if thermalization and 
hydrodynamic behavior set in later, they will be able to build only on a 
significantly reduced spatial eccentricity e x , and the resulting elliptic flow 
response will be correspondingly smaller. To reach a certain measured value 
of i>2 at a given impact parameter thus requires thermalization to set in be- 
fore free radial motion has reduced the spatial deformation so much that 
even perfect hydrodynamic motion can no longer produce the measured 
momentum anisotropy. This consideration yields a rigorous upper limit for 
the thermalization time r oqu . 

The dilution of the spat ial e ccentricity by collisionless radial free- 
streaming is easily estimated j ^ l ^ l using the analytic solution of the col- 
lisionless Boltzmann equation for the distribution function f(r,pr,r) of 
initially produced approximately massless partons (we only consider their 
transverse motion): 

/(r, p T , t+At) = f(r- cAt e p , p T ,r) . (24) 

Here e p is a unit vector in direction of pr- With Eq. H24JI it is straightfor- 
ward to compute the time-dependence of the spatial eccentricity: 

e (to+At) = / dx d y(y 2 ~ x2 ) I d2 PT f(r- cAt e p ,p T , r ) 
° / dxdy(y 2 +x 2 ) J d 2 p T f(r - cAt e p ,p T ,r ) 

J dxdy pTdpTdcpp [(y+cAr sin ip p ) 2 — (x+cAt cos <p P ) 2 ] f(r,pr, To) 
J dxdyp T dpTdip p [(y+cAr sin ip p ) 2 + (x+cAt cos tp p ) 2 } /(r, p T , r ) ' 

The initial distribution at To is even in x and y, and the initial transverse 
momentum distribution can be assumed to be locally isotropic. From this 
it follows directly that 

£2(70+ At) 



o(t ) 



(cAt) 
14 V 



r 2 ) 



(26) 



where (f 2 ) r is the azimuthally averaged initial transverse radius squared 
of the reaction zone. Inserting typical values for, say, Au+Au collisions at 
b = 7 fm one finds that a delay of thermalization by At = 2.5 fm/ c (3.5 fm/c) 
leads to a decrease of the spatial eccentricity by 30% (50%), without gen- 
erating any momentum anisotropy. The elliptic flow signal resulting from 
subsequent hydrodynamic expansion would then be degraded by a similar 
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percentage. Since at 6 = 7fm the RHIC data exhaust the hydrodynamic 
limit calculated with the full initial eccentricity e x (ro) at least at the 80% 
level, the thermalization time T eqn cannot be larger than about 1.75 fm/c. 

Excitation function of elliptic Row: 

At RHIC energies, the spectator nucleons (i.e. those nucleons in the two 
colliding nuclei which do not participate in the reaction) leave the reaction 
zone at midrapidity before the transverse dynamics begins to develop. This 
is why the spectator matter could be completely ignored in the initialization 
of the thermodynamic fields (see Sec. 12.4(1 . At lower energies, the spectator 
matter plays an active role in the dynamics as it blocks the transverse flow 
of mat ter into the reaction plane. Instead, its pressure squeezes the matter 
ou perpendicular to the reaction plane (i.e. in y-direction), resulting in 
a negative elliptic flow signal. A transition from negative to positive elliptic 
flow has been observed at the AGS in Au+Au collisions at beam ener- 
gies of 4-6 GeV_ per nucleonEHHI 

Relativistic hadron transport model 
calculation :) " I ' I indicate a need for a softening of the equation of state 
in order to quantitatively reproduce the data. 

Here we focus on the excitation function of the in-plane elliptic flow 
from the high end of the AGS energy range to RHIC and beyond. The left 
panel of Fig. 1241 shows the pr-hrtegrated elliptic flow for Pb+Pb or Au+Au 
collisions at fixed impact parameter b = 7 fm from a purely hydrodynamic 
calculation^ The excitation function is plotted versus the final particle 
multiplicity since hydrodynamics provides a unique relation between vi 
and dN/dy at fixed impact parameter but cannot predict the dependence 
of the latter on the collision energy. The horizontal arrows indicate ex- 
pected multiplicity ranges for RHIC and LHC before the first measurements 
were performed. We now know that Au+Au collisions at y / s = 200AGeV 
and 6 = yield charged particle multiplicity densities dN ch /dy~650 at 
midrapidityj ^ ^ somewhat lower than most people originally expected! ^'^ ) 
It is an easy exercise to translate this into total hadron multiplicity density 
at b = 7 fm as used in the Figure (roughly you have to multiply the above 
number by 3/4). El 

A characteristic feature of the hydrodynamic excitation function of V2 
(solid line in Fig. I24|l is a pronounced maximum around multiplicities corre- 
sponding to SPS energies ! ^ ^ | foll owed by a minimum in the RHIC energy 
domain (the 200AGeV datJEH 

correspond almost exactly to this mini- 
mum). The origin of this structure is the quark-hadron phase transition 
in the equation of stated It softens the matter in the phase transition re- 
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Pb+Pb, b=7fm 




total pion multiplicity density at y=0 

(GeV) 



Fig. 24. Left: Excitation function of the elliptic (solid) and radial (dashed) flow for 
Pb+Pb or Au+Au collisions at b = 7 fm from a hydrodynamic calculationBlThe collision 
energy is parametrized on the horizontal axis in terms of total particle multiplicity density 
dN/dy at this impact parameter. Right: A compilation of V2 data vs. collision e nergy 
from midcentral (12—34% of the total cross section) Pb+Pb and Au+Au collisionsli22J 

gion (small speed of sound c s ) which inhibits the buildup of transverse flow 
(both radial and elliptic). This effect is strongest for RHIC energies where 
a crucial part of the expansion when the source is still spatially deformed 
is spent in or close to the phase transition. At SPS energies most of the 
elliptic flow is generated in the hadronic phase (where w i), whereas at 
LHC energies essentially all elliptic flow is generated in the QGP phase 

Unfortunately, the available dataj I shown in the right panel of Fig. 1231 
do not support this hydrodynamic prediction. Given the success of the hy- 
drodynamic approach at RHIC energies, as pointed out in the earlier parts 
of this Section, this suggests a breakdown of hydrodynamics at lower colli- 
sion energies. One can think of at least two reasons for such a breakdown: 

(i) Lack of early thermalization: The equilibration times at SPS energies 
might be larger than assumed in the calculation (r cqu = 0.8 fm/c), so the 
hydrodynamic evolution starts later and builds on an already degraded spa- 
tial eccentricity. Given the smaller particle densities and collision rates at 
lower collision energies this possibility requires serious consideration. 

(ii) Lack of late thermalization: In the hydrodynamic calculations for the 
SPS, the spatial eccentricity does not disappear until significantly after the 
fireball matter has been fully converted to hadrons. More than half of the 
finally observed elliptic flow is thus generated during the hadronic staged If 
during this stage the hydrodynamic evolution is replaced by a microscopic 
kinetic description, the hadronic growth of the momentum anisotropy is sig- 
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nificantly reduced " ! I due to viscous effects in the hadronic cascade. In 
fact, such a combination of a hydrodynamic description until hadronization 
followed by a hadronic cascade afterwards leads to a monotonously increas- 
ing excitation function, in qualitative agreement with experimental How- 
ever, even in this hybrid approach the successful description of the SPS data 
still requires rapid thermalization to produce sufficient anisotropics already 
in the early collision stag this argues against possibility (i) above. 

The hybrid approach also does a better job than pure hydrodynamics in 
very peripheral collisions at RHIC energies where it produces a weaker re- 
sponse to the remaining spatial deformation during the hadronic stage than 
the pure hydrodynamic approach which overestimates i>2 at large impact 
parameters. 

If incomplete thermalization, especially somewhat later in the collision, 
is the reason for the failure of hydrodynamics at SPS energies, it would be of 
enormous help if we could perform elliptic flow studies with larger deformed 
systems than those created in peripheral Pb+Pb collisions. Fully central 
collisions between uranium nuclei offer such a possibility! 1 1 At similar 
spatial deformation as 6 = 7fm Pb+Pb collisions, the fireball formed in a 
side-on-side U+U collision is about twice as big in the transverse direction, 
a larger fraction of the elliptic flow is created before hadronization, and the 
fireball decouples several fm/c later and with smaller transverse velocity 
gradients^ All of these aspects should significantly improve the chances 
for a successful hydrodynamic description of elliptic flow in central U+U 
collisions down to SPS energies. This might open the door for confirming the 
hydrodynamic prediction of a non-monotonic behavior of vi as a function 
of collision energy which would unambiguously signal the softening of the 
equation of state near the hadronization phase transition^ 

Elliptic Row at non-zero rapidity: 

The hydrodynamic results presented so far were obtaine d with a code 
which explicitly implements longitudinal boost invarianc d^^^ lll ( see 
Sec. I2.6|l . By doing so one gives up all predictive power for the rapidity 
dependence of physical observables. Even if longitudinal boost invariance 
may be a reasonable approximation around midrapidity, it surely breaks 
down close to the longitudinal kinematic limit, i.e. near the projectile and 
target rapidities. 

In order to overcome these limitations, more elaborate 3+1 dimen- 
sional hydrodyn amic ca lculations have recently been done, mostly by two 
Japanese groups.^ 4 ! 141 ! 14 ^ 4 ^ 44 ! These require, of course, initial condi- 
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tions along the entire longitudinal axis. Since the re is relatively little hy- 
drodynamic evolution in the rapidity directional (due to the logarithmic 
nature of this variable), the initial longitudinal distributi on is rather tightly 
constrained by the measured final rapidity distributions (Mil xhe underly- 
ing assumption is, of course, local thermalization at all rapidities, which 
needs to be tested. 

Such a model can then predict the dependence of elliptic flow on rapid- 
ity. The PHOBOS collaboration has measured the elliptic flo w of uniden- 
tified charged particles over a wide range of pseudorapidityfSSI J n Fig- 
ure [23 these data are compared to a fully three-dimensional hydrodynamic 
calculation!^ The two different curves use different equations of state in 

0.08 



0.06 



^ 0.04 
0.02 



-6 -4 -2 2 4 6 

n 

Fig. 25. px-integrated elliptic flow for minimum bias Au+Au collisions at 
\/s = 130 A GeV as a function of pseudorapiditypHl compared with data from PHOBOS 
and 

the hadronic phase, either assuming chemical equilibrium (CE) until kinetic 
freeze-out or assuming decoupling of the abundances of stable hadrons al- 
ready at hadronization (PCE). One sees that as soon as one moves away 
from midrapidity the nice agreement of the data with the hydrodynamic 
predictions is lost. The origin of this is not yet fully clarified but is presum- 
ably a combination of both longitudinal boost invariance and local thermal 
equilibrium breaking down away from midrapidity. 

It may be worth commenting on some of the features of the data 
and the calculations: The slight peak at r] — (which is more strongly 
expressed in the data than in the calculation) is largely a kinematic ef- 
fect which arises wh en on e plots a flat rapidity distribution as a func- 
tion of pseudorapidityP^I Preliminary t>2 (rj) data from the 200 GeV run at 




February 5, 2008 7:22 



WSPC/Trim Size: 9in x 6in for Review Volume 



Hydrodynamic description of ultrarelativistic heavy-ion collisions 57 

appear to be less peaked in the region \rj\ < 1 and more consistent 
with the hydrodynamic calculations in this 77-domain. The hydrodynamic 
peaks at 77 ~ ±3 reflect mostly the initial conditions of the calculation, in 
particular the decreasing initial entropy density at forward and backward 
rapidities which, according to the discussion of the left panel of Fig. 1241 
should lessen the softening effects of the quark-hadron transition and in- 
crease the elliptic flow before dropping to zero at even larger rapidities. 
That this hydrodynamic feature is again not seen in the data may have 
similar reasons as the observation of a smaller than predicted elliptic flow 
at SPS energies (see previous subsection). 

The rapidity-dependent baryon and antibaryon spectr show that 
the net baryon number and thus the baryon chemical potential fib increases 
with rapidity \y\. With observables at finite rapidity one might therefore 
be able to explore the equation of state over a larger parameter space. 
This is of particular interest since recent lattice calculation s)!^ ! provide 
evidence for a tricritical point somewhere at fib > where the quark-hadron 
transition changes from a rapid crossover into a first-order phase transition. 
The effects of such a tricritical point in a hydrodynamically evolving system 
could lead to interesting signals and phenomena, such as bubble formation 
accompanied by large fluctuations! ^ 

Beyond hydrodynamics - elliptic How at high p^: 

Hydrodynamic flow flattens the transverse momentum spectra, and el- 
liptic flow flattens them more strongly in p x than in p y direction. As a 
result, hydrodynamics predicts v?. = (cos 2(p„) to grow monotonously with 
Pt , approaching V2 = 1 at px = 00 

In reality, high-px particles do not behave hydrodynamically since they 
escape from the fireball before having suffered a sufficient number of rescat- 
terings to thermalize their momenta. Whereas the hydrodynamic spec- 
tra drop exponentially at large pt, with an asymptotic slope reflecting a 
blueshifted freeze-out temperature^IHSI j" cff — j' dcc ^/l±^ (where v r here is 
the lar gest radial flow ve locity at any point inside the fireball), the measured 
spectrd ^W^l^^l^ l exhibit a power law tail at high px which increases 
the yield of high-px hadrons much above the hydrodynamic expectation. 

As reviewed elsewhere in this volume Qii^J these "hard" hadrons reflect 
production and medium modification mechanisms which, at sufficiently high 
Pt, can be computed directly from QCD. Above pr Z 4 GeV/c they domi- 
nate hydrodynamic particle emission ! I In the cross-over region between 
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the hydrodynamic domain (j?x$2GeV/c) and the perturbative QCD do- 
main (pt > 5 jzlOGsV/c) quark- coalescence may be the primary production 
mechanism! < | l^o | lU4 | lot> | lo ( \ j n ^ s pi c t ure; intermediate-px hadrons 
form by recombination of two or three partons with lower px and thus 
still partially reflect the hydrodynamic behavior which fully controls soft 

hadron productio n. 

Recent studie d 104 ' 156 ' suggest that at px > 5 GeV/c hadron production 
may start to become dominated by fragmentation of hard partons whose 
production can be calculated perturbatively. In this domain one again ex- 
pects momentum anisotropics which are correlated with the reaction plane, 
but from an entirel y diff erent effect which has nothing to do with collective 
hydrodynamic flowP^ High energy partons traveling through a medium 
with dcconfincd color charges are expected to lose energy via induced gluon 
radiation (see the accompanying article by Gyulassy et alP22 for details 
and references). In noncentral collisions, high momentum particles trav- 
eling along the short direction of the overlap zone (i.e. into the reaction 
plane) will lose less momentum and thus escape with higher px than par- 
tons emitted perpendicular to the reaction plane which have to cross a 
longer stretch of hot and dense matter. This leads to a positive v alue of 
V2 even at high px wich drops to zero logarithmically as px —> oo.' 1 ^ ' A 
qualitatively similar beha vior at high px is expected in the "Color Glass 
Condensate" pictureP^' 
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panel shows results from STAH v l from Au+Au collisions at 130 GeV, together with 
hydrodynamioEOand perturbative jet quenching calculations assuming various initial 
gluon densities. The right panel shows preliminary results of the elliptic flow of all charged 
hadrons from 200 A GeV Au+Au collisions for four centrality classes, as measured by the 
STAR collaborationlinSl 

Quantitative calculations of this effect, including the interpla y^ 1 ' with 
hydrodynamic elliptic flow at low px and the diluting effects of hydrody- 
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namic expansion on the matter traversed by the hard partoiff^EHH a t 
high pt, were recently performed. Generally, the observed signals (shown 
in Fig.l26|l are considerably larger than predicted. Even in the extreme limit 
of complete jet quenching (zero mean free path for hard partons), wher e V2 
reduces to an albedo effect whic h can be calculated geometricallypSI the 
predicted signal is still too smallP21 This suggests that in the px-region 
covered by Fig. 1261 additional contributions, other than parton energy loss, 
still contribute to the data. This may include unwanted two-particle corre- 
lations which are not correlated with the reacti on pl ane and which may be 
removable by a higher-order cumulant analysi J * once data with better 
statistics become available. 



4.2. Space-time information from momentum correlations 

Obviously, the fireballs formed in heavy-ion collisions are too small and 
shortlived to obtain a picture of their size and shape in coordinate space 
by usual methods, i.e. by scattering something of them and observing the 
diffraction pattern. Quantum statistical and final state interaction induced 
correlations among the momenta of the produced particles offer an alterna- 
tive approach to extract space-time inform ation about the emitting source. 
This tool, known as intensity interferometrypSluses the fact that the prob- 
ability of finding two particles with given pair and relative momenta in the 
same event is not simply the product of the independent probabilities to 
find each particle with the corresponding momenta, but reflects correla- 
tions between these momenta which are sensitive to the distance between 
the two particles when they were emitted. Although this tool doe s not allo w 
for a complete model-independent reconstruction of the sourceP ^^^ I it 
puts tight constraints on its spatial and temporal structure which, when 
combined with theoretical constraints for a consistent dynamical evolution 
and with supplementary experimental information from the single-particle 
spectra on the momentum-space str ucture of the source, allow to extract 
very detailed space-time information! 

An up-to-date review of this tool and a description of the underlying 
theoretical formalism can be found elsewhere in this volume in the article by 
Tomasik and Wiedemann!^ Here we only summarize the most important 
predictions for two-particle Bose-Einstein correlations among pions from 
the hydrodynamic model, in particular for non-central collisions where a 
new type of intensity interferometric analysis can complement the above 
discussion of momentum anisotropies generated during the early collision 
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stages with information about spatial deformations at the end of the colli- 
sion. Since any dynamical model relates these two aspects in a unique way, 
emission- angle dependent HBT-interferometry ( wher e the acronym stands 
for the initials of the originators of this methocP^) has the potential of 
putting much stronger constraints on these models than a study of the 
single-particle spectra or of central collisions alone. 

4.2.1. The hydrodynamic source function 

HBT interferometry is based on a fundamental relatiorPmaSHI (see 
Eq. H33fl below) which relates the two-particle correlation function to a 
Fourier transform of the source function (a.k.a. emission function) Si(x, K) 
for particles of species i. This emission function is a Wigncr density which, in 
hydrodynamic simulations, is replaced by its classical analogue, the phase- 
space probability density for finding a hadron i emitted from space-time 
point x with four-momentum K. If freeze-out is implemented in hydro- 
dynamics by sudden decoupling on a sharp freeze-out h ypersurfa ce £ as 
described in Sec. 12.51 the emission function takes the forrrP^^^OI 

S . (xK) 9i f K-d 3 a{x')5\x-x') 

W ' (2^ J^exp{[K-u(x')~^(x')]/T dec (x')}±l- [ ' } 

Phenomenological fits to spectra and HBT data often use a generalization of 
this form which replaces the ^-function by allowing for a spread of emission 
times ("fuzzy freeze-out" ) ) 171 | 172 | 173 | 

For a longitudinally boost invariant source (see Sec. 12. 6|) the freeze-out 
hypersurface can be parametrized in terms of the freeze-out eigentime as a 
function of the transverse coordinates, Tf(x, y). The normal vector d 3 a^ on 
such a surface is given by 

d 3 a— lcosh77, V±Tf(x,y), sinhr? J t/(x, y) dx dy di] . (28) 



With the four momentum — (Mr coshF, Kt, Mt sinhV), where Y 
and Mr are the rapidity and transverse mass associated with K, Eq. (|27() 
becomes 

/oo 
dr)dxdy[M T cosh(Y - V )-K T ■ V ± T f (x,y)] 

xf(K-u(x),x)S(r-T f (x,y)) (29) 

with the flow-boosted local equilibrium distribution /(if-u(a;), x) from H10JI . 

With this expression we can now study the emissivity of the source 
as a function of mass and momentum of the particles. For the purpose of 
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presentation we integrate the emission function over two of the four space- 
time coordinates and discuss the contours of equal emission density in the 
remaining two coordinates. We begin with calculations describing central 
Au+Au collisions at ^/snn = 130 GeV and focus on directly emitted pions, 
neglecting pions from unstable resonance decays. Resonance decay pions 
are known to produce non-Gaussian tails in the spatial emission distri- 
bution, increasing its width, but these tails are not efficiently picked up 
by a Gaussian fit t o the wi dth of the measured two-particle momentum 
correlation function .^^^23 A comparison of the experimental HBT size 
parameters extracted from such fits with the spatial widths of the emission 
function is thus best performed by plotting the latter without resonance 
decay contributions. 





K ± = 0.5 GeV 




5 10 
x (fm) 



Fig. 27. Pion source function S(x, K) for central Au+Au collisions at y 7 snn = 130 GeV. 
The upper row shows the source after integrating out the longitudinal and temporal coor- 
dinate, in the lower row the source is integrated over the longitudinal and one transverse 
coordinate (y). In the left column we investigate the case K = 0, in the right colu mn t he 
pions have rapidity Y = and transverse momentum Kt = 0.5 GeV in x directional 



Figure 1771 shows equal density contours at 10, 20, ... , 90% of the max- 
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imum in a transverse cut integrated over time and 77 (top row) and as a 
function of radius and time integrated over rj and the second transverse 
coordinate (bottom row). The dashed circle in the top row indicates the 
largest freeze-out radius reached during the expansion, the dashed line in 
the bottom row gives the freeze-out surface Tf(x,y=0) = tf(x,y—z—0). 
Pions with vanishing transverse momentum (left column) are seen to come 
from a broad region symmetric around the center and are emitted rather 
late. Pions with Kt = 0.5GeV pointing in x-direction, on the other hand, 
are emitted on average somewhat earlier and only from a rather thin, cres- 
cent shaped sliver along the surface of the fireball at its point of largest 
transverse extension. The reason for this apparent "opacity" (surface dom- 
inated emission) of high-px particles is that they profit most from the radial 
collective flow which is largest near the fireball surface. Low-px pions don't 
need the collective flow boost and are preferably emitted from smaller radii 
(where the flow velocity is smaller) when the freeze-out surface eventually 
reaches these points during the final stage of the decoupling process. 
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Fig. 28. Longitudinal cuts through the pion emission function S(x, K), integrated over 
transverse coordinates, for Y = pions with two different values of as indicated 1221 



Fig. [2H1 displays the time structure of the emission process along the 
beam direction. Especially for low tansverse momenta one clearly sees a very 
long emission duration, as measured in the laboratory (center-of-momentum 
frame). The reason is that, according to the assumed longitudinal boost in- 
variance, freeze-out happens at constant proper time r = t 2 — z 1 and extends 
over a significant range in longitudinal position z. This range is controlled 
by the competition between the longitudinal expansion velocity gradient 
(which makes emission of Y = pions from points at large z values un- 
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likely) and the thermal velocity smearing. If freeze-out happens late (large 
t), the longitudinal velocity gradient ~ 1/r is small and pions with zero 
longitudinal momentum are emitted with significant probability even from 
large values of \z\, i.e. very late in coordinate time t. Note that this is also 
visible in the lower left panel of Fig. 1271 where significant particle emission 
still happens at times where the matter at z = has already fully decou- 
pled. We will see shortly that this poses a problem when compared with 
the data. The long tails at large values of \z\ and t can only be avoided 
by reducing t/ (thereby increasing the longitudinal velocity gradient and 
reducing the z-range which contributes Y = pions) and/or by additionally 
breaking longitudinal boost-invariance by reducing the particl e de nsity or 
postulating earlier freeze-out at larger space-time rapidities \rj\ PH 

The spatial and temporal extensions of the emission process are char- 
acterized by the "spatial correlation tensor" 

S^(Y,K T ,^;b) = (x^x u ) (30) 

where i M (x^) is the distance to the center (x) of the emission region 

for momentum K and $ = Z(K^, b) is the azimuthal emission angle relative 
to the reaction plane. The averages are taken with the source function, 

/ s I d 4 xg(x,y,z,t)S(x,K) 

( 9 (x, y, z, t)) (K) = S # xS{XtK) • ( 31 ) 

The components of the spatial correlation tensor quantify the emission re- 
gions in terms of their spatial and temporal widths. These are directly 
related to the HBT size parameters extracted fr om t he width of the two- 
particle correlation function in momentum spaceP"3 

4.2.2. HBT-radii - central collisions 

Intensity interferometry is based on the analysis of the two-particle mo- 
mentum correlation function 

dN 



%^) = ^fi' (32) 



Quantum statistical effects (wave function (anti)symmetrization) between 
identical particles and final state interaction corrections between identi- 
cal or non-identical particles cause this correlation function to deviate 
from unity for small momentum differences P\~P2- One therefore con- 
veniently expresses it in terms of the momentum difference q — P\—p2 and 
the average momentum K = (pi+p 2 )/2. When these are supplemented by 
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the energy difference q° — E±— E^ and average energy K° = (E1+E2) /2 of 
the two particles, the resulting four-vectors satisfy the orthogonality re- 
lation KnQ ^ = (jrjg— rn 2 ) /2 ( = for identical particles). Under standard 
assumption; ] ^ l 1 "" ! (such as the absence of final state interactions) the 
2-pion correlation function can be related to the pion emission function 
S(x,K): 



C{q,K)*l 



J d 4 xS(x, K)e l "- X 2 



J d i xS(x, K) 



(33) 



For central collisions this does not depend on the azimuthal emission an- 
gle $. In a Cartesian coordinate system where the out-, side- and long- 
directions are defined parallel to K^, perpendicular to K^, and in beam 
direction, respectively, the source formed in a central collision is reflection 
symmetric under x s — > — x s . Exploiting longitudinal boost invariance of the 
source by selecting a frame for the analysis which moves with the longitudi- 
nal pair velocity /3l = Kl/K° (Longitudinal Co-Moving System, LCMS), 
one finds that in Gaussian approximatio n 1 the correlation function can 
be completely characterized in terms of three HBT-radii which depend only 
on the magnitude of K^: 

C(q, K)nl + exp [-R 2 (K T )q 2 - R 2 (K T )q 2 - Rf(K T )qf] . (34) 

(Without boost invariance the exponent contains in general a fourth 
ternJIEil and all HBT radii depend additionally on the pair rapidity Y.) 
These HBT-radii are directly related to the following combination of com- 
ponents of the spatial correlation tensor S^ v : 

R 2 S {K) = (x 2 s ) , (35) 
R 2 a (K) = (x 2 ) - 2/3 T ( Xo i) + f3 2 (t 2 ) , (36) 
R 2 (K) = (z 2 ), (37) 

where /3t = Kt I K° is the transverse pair velocity. 

From earlier hydrodynamic calculation J 1 77 1 was expected that a fire- 
ball evolving through the quark-hadron phase transition would emit pions 
over a long time period, resulting in a large contribution /3 T (t 2 ) to the 
outward HBT radius and a large ratio R a /R s . This should be a clear signal 
of the time-delay induced by the phase transition. It was therefore a big 
surprise when the first RHIC HBT dat a 1 1 ^ l 1 7» l 

yielded R D /R S w 1 in the 
entire accessible Kt region (up to 0.7GeV/c). In the mea ntime this finding 
has been shown to hold true out to K T > 1.2 GeV/cffSOl 
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Fig. 29. H BT r adii from hyd rody namic calculation^^ (solid lines) together with data 
from STAeBZHI and PHENIXM^I The dotted lines give hydrodynamic radii calculated 
directly after hadronization whereas the other two lines refer to different assumptions 
about the initial conditions (see text). 



Figure l2*§l shows a comparison of the experimental data with results from 
hydrodynamic calculations P2I Clearly, a purely hydrodynamic description 
with default initial conditions (solid lines in Fig. I29|) fails to describe the 
measured HBT radii. The longitudinal and outward radii Ri and R are 
too large, R s is too small, and the T^x-dependence of both R a and R s 
are too weak in the model. Since for longitudinally boost-invariant sources 
Ri is entirely controlled by the longitudinal velocity gradient at freeze- 
out which decreases as 1/t/, making Ri smaller withi n th e hydrodynamic 
approach requir es le tting the fireball decouple earlier^ or breaking the 
boost invarianceP^I As noted above, a smaller Ri would also reduce the 
emission duration, i.e. (P), and thus help to bring down R Q , especially if 
freeze-out at nonzero rapidities \t)\ > happens earlier than at midrapidity. 

The fireball can be forced to freeze out earlier by changing the freeze-out 
condition (e.g. by imposing freeze-out directly at hadronization, see dotted 
line in Fig. I29JI . This generates, however, serious conflicts with the single- 
Alternatively, one can allow transverse 



particle spectra (see Sec. 14.1. If) 
flow to build up sooner, either by seeding it with a non-z ero value already 
at r oqu (short dashed line labelled "hydro with FS'EEEl) or by letting the 
hydrodynamic stage begin even earlier (e.g. at Tf or m = 0.2 fm/c, long-dashed 
line). The last two options produce similar results, but do not fully resolve 
the problems with the magnitudes of Ri and R Q . 

The situation may improve b y ta king also the breaking of longitudi- 
nal boost-invariance into account but a fully consistent hydrodynamic 
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description has not yet been found. In particular, the sideward radius R s 
is still too small and the i^T-depenences of both R s and R are still to 
weak.EEni . g ■j agx ^ ^ Q gee Yxow to increase R s without also increasing R a 
and Ri which are already too large. Hybrid calculations^! in which the hy- 
drodynamic Cooper-Frye freeze-out is replaced by transition to a hadronic 
cascade at T crit , followed by self-consistent kinetic freeze-out, tend to in- 
crease R s by making freeze-out more "fuzzy", but at the expense of also 
increasing R Q and Ri in a disproportionate manner, mainly due to an in- 
crease in the emission duration. This makes the problems with the R /R s 
ratio even worse. 

In the past, a stro ng if j-d ependence of R s has been associated with 
strong transverse flowPSSQSZI It is therefore surprising that even the hy- 
drodynamic model with its strong radial flow cannot reproduce the strong 
ifx-dependence of R s measured at RHIC. Also, according to Eq. I|35|l the 
difference between R 2 and R 2 S can be reduced, especially at large i^x P^ 
if the positive contribution from the emission duration (t 2 ^ is compen- 
sated b y "source op acity" , i.e. by a strongly surface-dominated emission 
process! ^ 1 ^ 1 ^ In this case the geometric contribution (x 2 ) to R 2 is 
much smaller than R 2 = (x 2 .). Again, the source produced by the hydro- 
dynamic model (se e top right panel in Fig. l2"?)l is about as "opaque" as 
one can imagine j-^ and it will be difficult to further increase the difference 

This leaves almost only one way out of the "HBT puzzle" , namely the 
space-time correlation term — 2/3t {x a t) in expression Q35|l for R 2 . It corre- 
lates the freeze-out position along the outward direction with the freeze-out 
time. The hydrodynamic model has the generic feature that, in the region 
where most particles are emitted (see Fig. I27JI . these two quantities are 
negatively correlated, because the freeze-out surface moves from the out- 
side towards the center rather than the other way around. Hence the term 
— 2/3t (x t) is positive and tends to make R 2 larger than R 2 . The small 
measured ratio R /R s < 1 may instead call for strong positive x —t correla- 
tions, implying that particles emitted from larger x values decouple later. 
Hydrodynamics can not produce such a positive x —t correlation (at least 
not at RHIC energies). On the other hand, there are i ndications that micro- 
scopic models, such as the AMPT22SI an d MPcPSl models, may produce 
them, for reasons which are not yet completely understood. 

One should also not forget that Fig. [23] really compares two different 
things: The data are extracted from the width of the 2-particle correlator in 
momentum space while in theory one calculates the same quantities from 
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the source width parameters in coordinate space. The two produce iden- 
tical results only for Gaussian sources. The hydrodynamic source function 
shown in Figs. E3 and [2H1 are not very good Gaussians and show a lot of 
additional structure. We have checked, however, by explicit computation 
of the momentum-space correlation function using Eq. (EH that the non- 
Gaussian effects are s mall. b The largest non-Gaussian effects are seen in the 
longitudinal radius -B; P^ but although the corresponding corrections go in 
the right direction by making the Ri extracted from the momentum-space 
correlator smaller, the effect is only a fraction of 1 fm and not large enough 
to bridge the discrepancy wi th the data . 

It was recently suggesteoPHHaSi that neglecting dissipative effects 
might be at the origin of the discrepancy between the purely hydrodynamic 
calculations and the data. A calculation o f first -order dissipative corrections 
to the spectra and HBT radii at freeze-out EH with a "reasonable" value for 
the viscosity, yielded a significant decrease of Ri along with a correspond- 
ing strong reduction of the emission duration contribution to R , both as 
desired by the data. There was no effect on the x — t correlations, however, 
and only a weak effect on R s which went in the wrong direction, making 
it even flatter as a function of Kt- The rather steep i^T-dependence of the 
data for both R s and R and the larger than p redicted size of R s at low Kt 
are therefore not explained by this mechanism ! - ^ 1 l EEl Further more , th e el- 
liptic flow i>2 has been shown to be very sensitive to viscosity! •^^^ and 
th e visc osity values needed to produce the desired reduction in Ri turned 
ouJ l^ l to reduce i>2 almost by a factor 2, incompatible with the data. The 
"RHIC HBT puzzle" thus still awaits its resolution. 



4.2.3. HBT with respect to the reaction plane 

The discussion in Sec. 14.11 showed that azimuthal asymmetries in the 
momentum-space structure helped significantly in eliminating ambiguities 
of the global dynamical picture extracted from central collisions alone. Sim- 
ilarly, we can expect tight additional constraints from an analysis of spa- 
tial anisotropies, using HBT interferometry as a function of the azimuthal 
emission angle relative to the reaction plane in non-central collisions. In 
particular, dynamical models provide a characteristic relation between the 
momentum anisotropies generated early in the collision and the left-over 



b They would have been larger if we had included resonance decay pions in the emission 
function, as found by Lin et al. ) -^7v> l which was our main reason for not doing so. 



February 5, 2008 7:22 



WSPC/Trim Size: 9in x 6in for Review Volume 



68 Hydrodynamic description of ultrarelativistic heavy-ion collisions 

spatial deformation of the source at the end of the collision. This relation 
can be tested experimentally by combining elliptic flow measurements with 
the azimuthal dependence of the two-particle correlator and the extracted 
HBT radii, and this test can be used to ascertain the validity of the as- 
sumptions underlying the specific dynamical model. As an example, hydro- 
dynamics makes clear predictions for the time-dependence of the transverse 
spatial and momentum anisotropics, and, if valid, it should give the cor- 
rect sign and magnitude of the spatial deformation at freeze-out for the 
same set of model parameters which successfully reproduces the radial and 
elliptic flow in momentum space. Specifically, comparing the final spatial 
deformation to the initial one puts a constraint on the time between the 
beginning of transverse expansion and freeze-out, and this in turn may shed 
light on the origin of the discrepancies between the predicted and measured 
Ri values in central collisions discussed in the previous subsection. 

The formalism for HBT interfe rometry rela tive to the reaction plane 
was developed by Wiedemann et al ffSSHSSHSDl arK } is reviewed elsewhere 
in this volume Due to the lack of azimuthal symmetry in non-central 
collisions, there is no x s — > —x s symmetry and the exponent in Eq. (|34() 
contains all six terms, J^i j= s i R-ijQiQj) where the "HBT size parameters" 
Rij(Y, Kt, $) now also depend on the azimuthal emission angle 4> between 
the transverse emission direction Kt and the impact parameter b. For lon- 
gitudinally boost-invariant sources the transverse-longitudinal cross terms 
R 2 sl and R? ol still vanish in the LCMS frame, but there is an important out- 
side cross term R? os which is related to the spatial deformatio n of the sour ce 
in the transverse plane. For 6^0 equations (|35|l generalize tc™HHHni 

R 2 S (®) = \{S XX + S yy )- - (S xx ~ Syy) cos2$ -S xy sin 2$ (38) 

= ^(S XX + Syy) + -{S XX - Syy) COS 2$ + S XJ/ SU1 2$ 

-2p T (S tx cos $ + S ty sin $) + PlSu (39) 
R 2 os ($) = S xy cos2$- hs xx - S' ra )sin2$ 

+/3 T (5 te sin$-5 ty cos$) (40) 

Rf($)=S zz . (41) 

Here the spatial correlation tensor S^ v — is specified in reaction- 

plane coordinates (a: points along the impact parameter and y is perpen- 
dicular to the reaction plane) where, due the reflection symmetry of the 
overlap zone with respect to the x and y axes, it is most easily evaluated. 
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The indicated explicit ^-dependence arises from the rotation between the 
outward direction £ ||_Kt and the x-axis. In addition, the components S^v 
of the spatial correlation tensor, being defined as expectation values with a 
^-dependent emission function S(x, K) = S(x, Y, if-r, contribute an im- 
plicit ^-dependence which is not shown here. Both types of ^-dependences 
can be analyzed together, exploiting the symmetries of the emission func- 
tion with resp ect to the reaction plane and to projectile-target interchange. 
One find J Ifhat in general R 2 , R 2 and R 2 are superpositions of cosines of 
even multiples of $ while R 2 S is a superposition of sines of even multiples 
of $. In lowest order R 2 , R 2 , and Rf thus oscillate with cos(2$) around 
some constant average while R 2 S oscillates with sin(2$) around zero. 




Fig. 30. Contours of constant emission density for the s ourc e created in 130j4GeV 
Au+Au collisions at b = 7 fm as predicted by hydrodynamics t^The top left panel shows 
the contours for = 0, the three other panels show contours for = 0.5 GeV pointing 
in 3 different directions. The thick solid oval lines indicate the maximum extension of 
the source. 



In SecEPwe saw for non-central collisions that, as time evolves, the ini- 
tial out-of-plane deformation of the nuclear reaction zone decreases, crosses 
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zero and eventually turns into an elongation along the impact parameter 
direction (see Fig.EJ. At different times the freeze-out surface (see Fig. l27|l 
thus reflects different spatial deformations. Hence, the effective deformation 
probed by the two-particle correlation function is the average of the spatial 
eccentricity taken over the freeze-out surface. FigureOOIshows contour plots 
of the emission function from the same hydrodynamic calculations which 
reproduced the spectra and elliptic flow data from 130^4 GeV Au+Au colli- 
sions at RHIC (see Sec. 14. If) , for an impact parameter b = 7fm. The under- 
laid outer contour shows that, at the time of its largest spatial extension, 
the hydrodynamic source is still elongated out-of-plane, although much less 
so than at the beginning. As seen in the upper left panel, this out-of-plane 
elongation is probed by low-i^x pions which, just as in the central colli- 
sion case shown in Fig. 1271 are emitted from the entire region inside this 
outer contour. In contrast, pions with non-zero transverse momentum are 
emitted from relatively thin slivers near the outer edge of the fireball and 
therefore do not directly probe the overall deformation of the source. The 
shapes of their emission regions are controlled by an interplay between the 
curvature of the outer edge of the source and the strength of the anisotropic 
transverse flow, which both vary with the azimuthal emission anglel^ 



JO 



-10 



IPES ^ 
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Fig. 31. Contours of constant emission density for a strongly in plane oriented sourceEH 
The left panel shows the contours for pions with = 0, the right panel shows contours 
for pions with Ki = 0.5 GeV emitted in 3 different azimuthal directions. In both panels 
the thick solid line indicates the maximum extension of the source. 



In order to see how these patterns would change if at freeze-out the 
source were actually elongated along the reaction plane, and how the HBT 
correlation function would reflect this difference, we show in Fig. 1311 the 
emission function for a hydrodynamic source which evolved from a much 
higher initial energy density. The calculation shown in Fig. |^is for Au+Au 
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collisions at b — 7im, assuming an initial central temperature of 2GeV at 
T e qu = 0.1 fm/c and decoupling at Tdcc = 100 MeVP^I With such initial con- 
ditions it takes the fireball much longer to reach freeze-out and, as seen in 
the Figure, it has enough time to develop a strong in-plane elongation before 
decoupling. (IPES stands for "in-plane elongated source".) The emission 
function for low-fCx pions is again seen to reflect the overall in-plane elon- 
gation of the source at freeze-out (left panel), whereas the emission func- 
tions for pions with transverse momentum = 500MeV/c (right panel) 
do not probe the overall source deformation, but rather the curvature of its 
outer edge as a function of the azimuthal emission direction. 

These pictures suggest that, if one wants to measure the overall spatial 
deformation of the source at freeze-out, one should concentrate on pions 
with small transverse pair momenta K^. This is confirmed by the plots 
in Fig. 1321 which show the azimuthal oscillations of the four non- vanishing 
HBT radius parameters for several values of K^. At Ki — the trans- 




40 


I 20 






\ 




20% 




c cc° 










200_ 


% 1 


1= 




100g 




"aT 


-1 




Fig. 32. Oscillations of the HBT radii for different transverse pair-momenta in RHIC 
collisions (left) and the source elongated into the reaction plane (right) EH The thin 
circled lines in the right panel show the geometric contributions to the HBT radius 
parameters. 



verse radius parameters R 2 ., R 2 and R 2 S show opposite oscillations for the 
out-of-plane (left) and in-plane (right) elongated sources, and the signs of 
these oscillations reflect the signs of the geometric source deformation as 
expected. For example, at RHIC energies R 2 S oscillates downward, implying 
a larger sideward radius when viewed from the x direction (i.e. within the 
reaction plane) than from the y-direction (i.e. perpendicular to the reaction 
plane). For the in-plane elongated source (IPES) some of the oscillation am- 
plitudes change sign at larger transverse momentum. This change of sign 
originates from an intricate interplay between geometric, dynamical and 
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temporal aspects of the source at freeze-out!^ 

At this moment the significance of the ifx-dependence of the oscillation 
amplitudes for the HBT radius parameters is still largely unexplored. In the 
hydrodynamic model it is crucially affected by the Cooper-Frye freeze-out 
criterium which strictly limits the source size and provides a sharp radial 
cutoff fo r the distribution of emission points. From the analysis of central 
collision Jl^ l it is known that the emission-angle averaged transverse HBT 
radii (in particular the outward radius R a ) exhibit stronger Kt dependence 
for sources with a sharp surface (such as a box-like density distribution) 
than for Gaussian sources. One might expect similarly strong differences for 
the i^x-dependence of the oscillation amplitudes in non-central collisions. 
Model calculations in which the surface diffuseness of the source can be 
varied confirm this expectation and show a particularly stron g effect of this 
parameter on the oscillations of R 2 Q as a function of jf T j l91 | 

The oscillation pattern shown in the left panel of Fig. 03 agrees 
qualitative ly with measurements by the STAR Collaboration at 130 and 
200 A GeVp22D231 Figure OS shows preliminary data from Au+ Au colli- 




sions at snn = 200 GeV for three different centrality classes! Although 
these data are integrated over and not yet final and thus should not 
be overinterpreted, one observes that the oscillation amplitudes follow the 
pattern in the left panel of Fig. |23 but not that in the right panel. This 
implies that the fir eball form ed in these collisions is still out-of-plane elon- 
gated at fxeeze-outpMnSl as predicted by hydrodynamics. Whether the 
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deformation extracted from the data is larger than that in the hydrody- 
namic model (which, if true, would point to earlier freeze-out as discussed 
at the end of the previous subsection) is not yet clear. Higher statistics data 
which are presently being analyzed will soon answer this question. 
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5. Summary and Conclusions 

In this review we have summarized some of the most recent theoretical and 
experimental results indicating that the fireball created in Au+Au collisions 
at RHIC thermalizes very quickly, evolves through an extended dynamic 
expansion stage which is governed by intense rescattering, follows the laws 
of ideal hydrodynamics, and probes the equation of state of nuclear matter 
at temperatures between 350 MeV and 100 MeV. The data are sensitive 
to the softening effects of a phase transition on the nuclear equation of 
state and were shown to be consistent with the existence of a quark-hadron 
phase transition at energy densities around lGeV/fm 3 . Equations of state 
without such a phase transition, such as an ideal gas of massless partons or a 
hadron resonance gas with approximately constant sound velocity w 1/6 
extending to very high energy densities, are experimentally excluded. 

The most important result of the comparison between theory and exper- 
iment presented in this review is very strong evidence that hydrodynamics 
works at RHIC, i.e. that it is able to describe the bulk of the data on soft 
hadron production at pt < 1.5 — 2GcV/c. This includes all hadronic mo- 
mentum spectra from central (b = 0) to semiperipheral (b < 10 fm) Au+Au 
collisions, including their anisotropics. Since the anisotropies are sensitive 
to the spatial deformation of the reaction zone and thus develop early in 
the collision, they provide a unique probe for the dense matter formed at 
the beginning of the collision and its thermalization time scale. The data 
on elliptic flow can only be understood if thermalization of the early par- 
tonic system takes less than about lfm/c. At this early time, the energy 
density in the reaction zone is about an order of magnitude larger than 
the critical value for quark deconfinement, leading to the conclusion that 
a well-developed, thermalized quark-gluon plasma is created in these colli- 
sions which, according to hydrodynamics, lives for about 5-7fm/c before it 
hadronizes. 

The mechanisms responsible for fast thermalization are still not fully 
understood. Perturbative rescattering of quarks and gluons was shown to be 
insufficient, leading to thermalization times which are a factor 5-10 too long. 
The quark-gluon plasma created at RHIC is thus strongly nonperturbative. 

As a result of the softening effects of the quark-hadron transition on 
the nuclear equation of state, hydrodynamics predicts a non-monotonous 
beam energy dependence of elliptic flow. In particular, the elliptic flow 
at SPS energies is expected to be larger than at RHIC, inspite of some- 
what smaller radial flow. This is not borne out by the existing data from 
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Pb+Pb collisions at the SPS, due to the unfortunate breakdown of the 
hydrodynamic model at the lower SPS energy. It is suggested to explore 
the excitation function of elliptic flow at RHIC by going down as much as 
possible with the beam energy, using central collisions between deformed 
uranium nuclei. These would, at similar spatial deformation, provide much 
larger reaction volumes than peripheral Au+Au collisions, improving the 
chances for hydrodynamics to work even at lower energies where the fire- 
balls are less dense and thermalization is harder to achieve. A decreasing 
elliptic flow as a function of increasing collision energy in U+U collisions 
at RHIC, followed by an increase as one further proceeds to LHC energies, 
would be an unmistakable signature for the existence of the quark-hadron 
phase transition. 

The gradual breakdown of hydrodynamics in peripheral Au+Au colli- 
sions between RHIC and SPS energies goes hand-in-hand with a gradual 
breakdown of hydrodynamics in Au+Au collisions at RHIC as one moves 
away from midrapidity. The elliptic flow was found to drop dramatically 
at pseudorapidities > 1, and hydrodynamic calculations arc unable to 
reproduce this feature. It is likely that both observations are related, and 
again central U+U collisions may clarify this issue. 

Whereas the momentum-space structure of the observed hadron spectra 
is described very well by the hydrodynamic model, at least for pt < 1.5 — 
2 GeV/c, this is not true for the two-particle momentum correlations which 
reflect the spatial structure of the fireball at hadronic freeze-out. The hy- 
drodynamic model shares with most other dynamical models a series of 
problems which have become known as the "RHIC HBT puzzle" . Typically, 
the longitudinal HBT radius is overpredicted while the sideward radius and 
its transverse momentum dependence are underpredicted. Furthermore, all 
calculations give an outward radius which is significantly larger than the 
sideward radius, in contradiction with experiment which shows that both 
radii arc about equal. 

This implies a serious lack of our understanding of the decoupling pro- 
cess in heavy-ion collisions and how its details affect our interpretation of 
the size parameters extracted from two-particle intcrfcrometry. However, 
these two-particle correlations are only fixed at the end of the collision 
when the hadrons cease to interact strongly with each other. Our failure to 
reproduce the correlations within the hydrodynamic model therefore does 
not affect our successful description of the elliptic flow patterns (which are 
established much earlier) within the same model. This is important since 
a quantitative understanding of the global collision dynamics, especially 
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during its early stages, is crucial for a quantitative interpretation of "hard 
probes" such as jet quenching, direct photon and dilepton production and 
charmonium suppression. These cannot be successfully exploited as "early 
collision signatures" without a proper understanding of the global collision 
dynamics which the hydrodynamic model provides. 

As shown in the last Section, this global picture can be further con- 
strained by emission-angle dependent HBT interferometry which allows to 
correlate momentum anisotropics measured by v-i with spatial anisotropics 
at freeze-out measured by the azimuthal oscillation amplitudes of the HBT 
radius parameters. Preliminary data confirm the hydrodynamic prediction 
that, at RHIC energies, the source is still slightly elongated perpendicular 
to the reaction plane when the hadrons finally decouple. At the LHC, the 
spatial deformation at freeze-out is expected to have the opposite sign, the 
fireball being wider in the reaction plane than perpendicular to it. 
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